Last edited: January 18, 2011 

Preprint typeset using I^T^jK style cmulatcapj v. 03/07/07 



UV-TO-FIR ANALYSIS OF SPITZER/IRAC SOURCES IN THE EXTENDED GROTH STRIP I: 
MULTI- WAVELENGTH PHOTOMETRY AND SPECTRAL ENERGY DISTRIBUTIONS 

Guillermo Barro 1 , P.G. Perez-Gonzalez 1,2 , J. Gallego 1 , M. L. N. Ashby 3 , M. Kajisawa 4 , S. Miyazaki 5 , V. 

VlLLAR 1 , T. YAMADA 4 , J. ZAMORANO 1 

(Dated: Submitted: January 18, 2011) 
Last edited: January 18, 2011 



ABSTRACT 

We present an IRAC 3.6+4.5 /j,m selected catalog in the Extended Groth Strip (EGS) containing 
photometry from the ultraviolet to the far-infrared and stellar parameters derived from the analysis 
of the multi-wavelength data. In this paper, we describe the method used to build coherent spectral 
energy distributions (SEDs) for all the sources. In a companion paper, we analyze those SEDs to 
obtain robust estimations of stellar parameters such as photometric redshifts, stellar masses, and star 
formation rates. The catalog comprises 76,936 sources with [3.6]<23.75 mag (85% completeness level 
of the IRAC survey in the EGS) over 0.48 deg 2 . For approximately 16% of this sample, we are able 
to deconvolve the IRAC data to obtain robust fluxes for the multiple counterparts found in ground- 
based optical images. Typically, the SEDs of the IRAC sources in our catalog count with more than 
15 photometric data points, spanning from the ultraviolet wavelengths probed by GALEX to the far- 
infrared observed by Spitzer, and going through ground- and space-based optical and near-infrared 
data taken with 2-8 meter class telescopes. Approximately 95% and 90% of all IRAC sources are 
detected in the deepest optical and near-infrared bands. These fractions are reduced to 85% and 70% 
for SNR>5 detections in each band. Only 10% of the sources in the catalog have optical spectroscopy 
and redshift estimations. Almost 20% and 2% of the sources are detected by MIPS at 24 and 70 /im, 
respectively. We also cross-correlate our catalog with public X-ray and radio catalogs. Finally, we 
present the Rainbow Navigator public web-interface utility, designed to browse all the data products 
resulting from this work, including images, spectra, photometry, and stellar parameters. 
Subject headings: galaxies: starburst — galaxies: photometry — galaxies: high-redshift — infrared: 
galaxies. 



1. INTRODUCTION 

Multi-wavelength observations of blank fields provide 
a fertile ground for studies of the evolution of galaxies 
from the early Universe. Indeed, in the past decade 
we have advanced amazingly in our knowledge about 
the formation of galaxies thanks to deep held imaging 
and spectroscopic surveys. The extraordinary success 
of these surveys is sustained by the coordinated effort 
of several telescope facilities, institutions, and research 
groups that gather large collections of multi-wavelength 
photometry and spectroscopy, providing the entire scien- 
tific community with a vast pool of data to analyze. Re- 
markable examples of this kind of projects are the Hubbl e 
Deep Field (HDF) observations ([Williams et al.lll996l ). 
the Classifyin g Objects by M edium- Band Observa- 
tions project (jWolf et all 120011 C OMBO 17), Great Ob- 
servatories Optical Deep Survey dGiavalisco et alJl2004L 
GOO DS), VIMOS-VLT Deep Survey dLe Fevre et al.1 
120051 VVDS), All- wavelength Extend ed Groth strip In- 
ternational Survey (jDavis et all 120071, AE GIS), or Cos- 
mic Evolution Survey (|Scoville et alll2007l COSMOS). 
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Nevertheless, the full scientific exploitation of these 
surveys unavoidably needs a consistent merging of the 
data coming from heterogeneous sources (with differ- 
ent depths and resolutions) to build catalogs of galaxies 
characterized with panchromatic photometry and spec- 
troscopy. Although a substantial effort has been devoted 
to create homogeneously processed multi- wavelength cat- 
alogs for the most important helds (e.g., Chandra Deep 
Field South, JWolf et al.ll2001l: Hubble D eep Fi eld North 
[ Steidel et al.l I2003L iCaoak et al.l l2004l: EGS. iCoil et"aL1 
12004 lllbert et al.l 127)06^ 7 many of these catalogs are 
selected in the optical bands (i.e., rest- frame ultra- 
violet -UV- at high redshift), and lack near- infrared 
(NIR) imaging (an important dataset for the stud- 
ies of galaxy populations at high redshift) reaching 
depths that match the optical observations. Fortu- 
nately, the proliferation of deep and wide NIR sur- 
veys has supported the publication of an increasing 
number of multi-ban d samples selected i n the if -band 
dGrazian et al.l 120061 iQuadri et ail 120071; IWuvts et all 
120081) or iRAC-ban ds dRowan-Robinson et all 120081 : 
iPerez-Gonzalez et all [2005. 200|, hereafter, PG05 and 
PG08). A detailed UV-to-NIR coverage of the spectral 
energy distribution (SED) improves the estimates of im- 
portant stellar parameters, such as t he mass, age, extinc- 
tion, or star form ation rate (see, e.g. JWalcher et ai1l2008l : 
ISalim et al.ll2007l) . Moreover, in order to obtain the most 
reliable stellar mass estimates for z^2-3 galaxies, and to 
distinguish young dusty starburst from quiescent galax- 
ies at these redshifts, we need to obtain data probing the 
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rest-frame NIR for t h ese populations (iBundv et al.ll2006l; 
Pozzetti et "all l2007t IWuvts et al.l 120071 : IWilliams et al.l 
2009tlMuzzin et al.ll2009t) . This became possible with the 
launch of Spitzer and th e use of one of its instruments, 
IRAC ([Fazio et al.ll2004| ). which covers the 3.6-to-8.0 ^m 
spectral range. 

The power of multi-band catalogs is significantly en- 
hanced when optical and NIR data are complemented 
with mid-IR (>5/im) to radio fluxes, such as those from 
Spitzer/MIPS, SCUBA, VLA or Herschel surveys. These 
data directl y probe the emission of the dust comp onent 
of galaxies fSantini et afll2009t IWuvts et al.l 12008). Ac- 
tively star-forming galaxies harbor large amounts of dust, 
which cause that a fraction of the UV emission, directly 
related to the ongoing star formation, is extincted and 
re-emitted in the IR. Thus, modeling the IR emission 
offers not only a complementary approach to estimate 
the star formation rate (SFR) of a galaxy but also an 
improved measurement of the intrinsic UV extinction 
(iReddv et al.ll2006b iDaddi et al.l l2007HSalim et~alll2007t 
llglesias-Paramo et alj 12007). The IR approach to stud- 
ies of the star formation becomes particularly relevant 
at higher redshifts, where the number of luminous in- 
frared galaxies (LIRGs, whose integrated IR luminosity 
is L(IR)>10 11 Lq) and their contribution to the cosmic 
SFR d ensity increase sign ificantly (Chary&Elbaz 2001, 
PG05. ICaputi et al.ll200l . 

The downside to having an exceptional data pool avail- 
able in many regions of the sky is that the data quality is 
largely heterogeneous. Unfortunately, the high-redshift 
community still lacks the existence of a unified database 
that facilitates the access to the multiple datasets and 
resources, similar to local extragalactic databases such 
as the NASA Extragalactic Databa se (NED) or the 
Sloan Digitalized Sky Survey (SDSS; lYork et al.l 1200(1 
database. 

The purpose of this work is to present a NIR selected 
sample of galaxies with well sampled spectral energy dis- 
tributions, and analyze their properties maximally ben- 
efiting from the panchromatic data. To do this, we 
have built an IRAC-3.6+4.5/im selected photometric and 
spectroscopic catalog including data from X-ray to ra- 
dio wavelengths for 76,936 galaxies at 0<z<4 in the Ex- 
tended Groth Strip (EGS). 

The EGS has been intensively observed as a part of 
the All- Wavelen gth Extended Gro th Strip International 
Survey (AEGIS : lDavis et al.ll2007t ) collaboration in order 
to assemble an exceptional multi- wavelength dataset, in- 
cluding deep optical imaging from the CFHTLS, HST 
coverage in two bands, UV data from GALEX, and mid- 
IR and far-IR photometry from Spitzer. In addition, 
the EGS is the key field' for the DEEP2 survey, one 
the largest and d eepest spectroscopic surveys to date 
([Davis et al.ll2003T ). with more than 10,000 optical spec- 
tra down to R~24. This vast dataset converts the EGS in 
one of the main fields for the study of galaxy evolution at 
different epochs of the lifetime of the Universe. Our goal 
is get advantage of this impressive multi- wavelength data 
collection and use it to build UV-to-FIR SEDs, whose 
analysis will allow us to obtain estimations of interest- 
ing parameters, such as the photometric redshifts, stellar 
masses, and SFRs. 

The photometric catalog along with the photometric 
redshifts and the inferred stellar parameters are intended 



to become a multi-purpose resource useful for many dif- 
ferent scientific goals. Some of them will be presented in 
forthcoming papers. We make all the catalogs publicly 
available through our website and through a dedicated a 
web-interface, dubbed Rainbow Navigator, conceived to 
facilitate the access to the data, but also to serve as a 
permanent repository for updates in these catalogs, or 
similar catalogs in other cosmological fields (e.g., those 
presented in PG08 for the GOODS fields). 

In this paper, we concentrate on the description of 
the dataset and the methods developed to measure the 
merged photometry for the IRAC sample in the EGS. 
We also analyze the multi-band properties of the sample 
to understand the main properties of the Spitzer sur- 
veys. All the multi- wavelength photometry is released in 
a public database, conceived to allow the astronomical 
community to access all the results from our work and 
use them for their own purposes. In a companion paper 
(hereafter Paper II) , we will present our methodology to 
fit the SEDs presented in this paper, and to estimate 
photometric redshifts, stellar masses, and SFRs out of 
them. We will also assess the quality of the inferred pa- 
rameters, analyzing in detail their intrinsic systematic 
and random uncertainties. 

The outline of this paper follows. In §[2] we present the 
available datasets that we have compiled for this paper. 
In §[3] we present the techniques developed to extract the 
IRAC 3.6+4.5 ^m selected sample and we discuss the 
properties of the IRAC photometry of the sample. In §|4] 
we present the methods developed to build the merged 
multi-band photometric catalog. In § [5] we describe in 
detail the photometric properties and reliability of the 
catalog. In § [6] we describe the format of the published 
catalogs, the database built to allow an easy access and 
handling of those, baptized as the Rainbow Cosmolog- 
ical Surveys Database, and the publicly available web- 
interface to surf the database, Rainbow Navigator. 

Throughout this paper we use AB magnitudes. We 
adopt the cosmology Hq = 70 km -1 s -1 Mpc -1 , Q m = 0.3 
and Q\ = 0.7. 

2. DATA DESCRIPTION 

The Extended Groth Strip (EGS; a = U h l7 m , 
5 = +52°30') is one of the most targeted cosmological 
deep fields. Noticeably, a comprehensive panchromatic 
dataset has be en compiled in thi s field within the AEGIS 
collaboration (|Davis et al.ll2007D . 

The sample of galaxies studied in this paper is based 
on an IRAC-3.6+4.5/j,to selection. This choice obeys to 
several reasons. First, the IRAC bands are specially tai- 
lored to probe the rest-frame near-infrared (NIR) fluxes 
of distant galaxies, thus being the per fect tool for stud- 
ies of massive galaxies at high- z (e.g.. iRodighiero et al.1 
I2007t iMancini et al"1l2009D . Second, the quality, depth 
and ubiquity of the IRAC observations in the so-called 
cosmological fields favors the assembly of coherent flux 
limited catalogs over large cosmological volumes (e.g., 
PG08), which is the cornerstone of the observational cos- 
mology. The IRAC bands offer an alternative to the less 
efficient NIR ground-based surveys, and provide a start- 
ing point to consistently anchor surveys at longer wave- 
lengths (Spitzer/MIPS, Herschel, A LMA). Despite its 
lower spatial resolution (FWHM~2"; iFazio et at] 12004(1 
compared to optical/NIR ground-based surveys, the im- 
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age quality of IRAC is very stable and several authors 
have been able to perform deblending techniques suc- 
cessfully mergin g IRAC catalogs into their p anchromatic 
datas ets (PG05. iGrazian et al.ll200l PG08, lllbert et al.l 
120091) . 

In the rest of this section, we describe the multi- 
wavelength datasets in EGS compiled for this paper to 
characterize the IRAC selected sample. Table Q] summa- 
rizes the main characteristics of these datasets, includ- 
ing the depth, the area and the image quality. Figure [T] 
shows the footprints of the various surveys, highlighting 
the area with higher band coverage. 

2.1. X-ray data 

The EGS region has been observed in the X-ray (0.5-10 
keV) using Chandra/ ACIS during two observation cycles 
in 2002 and 2005. The AEGIS-X survey covers an area 
of 0.67 deg 2 in 8 pointings (~17'xl7' each), completely 
overlapping with the region covered by IRAC. The nom- 
inal exposure time of the frames is 200 ks per pixel, 
reaching limiting fluxes of 5.3xl0~ 17 and 3.8xl0~ 16 
erg cm 2 s-1 in the soft (0.5-2 keV) and hard (2-10 keV) 
bands, respectively. In this work, we use the data reduc- 
tion and point source catalogs pub lished in lLaird et al.1 
(2009, see also iNandra et al.1 120051 for a first version of 
the catalogs). The two-band (soft and hard) merged cat- 
alog comprises 1325 sources with <1.5% spurious detec- 
tions. The authors identified optical and NIR counter- 
parts for 1013 and 830 sources, respectively, from the 
CFHT legacy survey optical catalog (see j j |2.3.Tj) . and 
the Spitzer /IRAC catalog of iBarmbv et~alT (|2Q08j). The 
cross-match is based on maximum likelihood method 
with a search radius of 2" (more than a factor of 1.5 
the rms of their astrometric accuracy) . 

2.2. Ultraviolet data 

The Ga l axy Evolution Explorer (GALEX, 
iMartin et all I2005D observed the EGS over three 
consecutive years providing deep UV data in two 
channels (FUV at 153 nm, and NUV at 231 nm) 
on a 1.13 deg 2 circular area around a = 14' l 20 m , 
S = +52°47'. The total exposure time for the composite 
stacks are 58 ks and 120 ks in the FUV and NUV filters, 
respectively. The approximate limiting magnitude in 
both bands is ~25.1 mag. 

2.3. Optical data 
2.3.1. CFHTLS/CFHT12K 

A 1° x 1° square region centered on the coordinates 
a = 14 ft 19 m 27 s , 5 = +52°40'56" was observed within the 
Canada-France Hawaii Telescope (CFHT) Legacy Survey 
(CFHTLS; sector D3). Using the MEGACAM camera on 
the CFHT, imaging data were obtained in five broadband 
filters over the wavelength range 350 nm<A<940 nm: 
u*, g' , r' , i', z'. The overall exposure times ranges 
from 4h to 40h (for the u* and z' bands, respectively), 
reaching limiting magnitudes between 26-27 mag. The 
reduced images have gone through several releases. The 
data used in this paper are part of the CFHTLS 'T0004' 
release produced at the TERAPIX data center (Gwyn 
et al. 2010, in preparation 2 ). Although the total mo- 
saicked area covers 1 deg 2 , the overlap with the IRAC 



observations is just 0.35 deg 2 . 

In addition to the CFHTLS data, shallower images 
were acquired in the B, R and I bands over a larger 
area, using the wide-field 12K mosaic camera on CFHT. 
These observations were intended to provide optical cov- 
erage for the DEEP2 spectroscopic survey, extending 
through the whole EGS in 4 different pointings. The ob- 
serving times were 3-6h, to achieve a limiting magnitude 
of ^25 mag. The photo metric catalogs w ere published 
within the DEEP2 DR1 (jCoil et al.ll2~00l . and the raw 
(non-reduced) images can be retrieved from the CADC 
archive 3 . We downloaded the raw and calibrations files 
from the archive to perform our own data reduction and 
cataloging. This reduction was carried out with the 
IRAF task mscred. The photometric calibration was 
performed using the BRI public catalogs. The average 
dispersion of the photometry comparison is smaller than 
0.03 mag. 

2.3.2. MMT/Megacam 

The EGS was ob served with the Megacam camera 
(|McLeod et al.l 12006ft on the MMT in the u'giz bands 
during several campaigns from 2006 to 2009. Four point- 
ings of the 24' x 24'Megacam field of view were arranged 
to cover the full extension of the 2 degrees long IRAC 
mosaic. 

The MMT/Megacam imaging data were reduced with 
a combination of standard IRAF/mscred routines and 
custom software. Final mosaics were calibrated against 
the coextensive SDSS photometry with appropriate color 
corrections. The 5<r limiting magnitudes within a 2" di- 
ameter aperture for point sources varies among the four 
different pointings but are approximately 26.5, 27.2, 26.0, 
and 26.0 mag in the u',g,i, and z bands, respectively. 
The effective seeing in the final mosaics is roughly 1" 
FWHM or slightly better in some bands. 

2.3.3. HST/ACS 

As a part of a GO program (PI: Davis), the HST Ad- 
vanced Camera for Surveys (ACS) acquired deep imaging 
of EGS in two optical bands: i^606IV and F8UW (here- 
after the Vqo6 and igi4 images). The ACS survey covers 
an area of ~710.9 arcmin 2 (10.1'x70.5') on 63 contiguous 
tiles following the direction of the IRAC mosaic, covering 
approximately 50% of the total area surveyed by IRAC, 
and ~ 80% of the area overlapping with the CFHTLS 
observations. The fully reduced, drizzled frames and 
calibration produ cts were released by the AEGIS Team 
(jLotz et al.M2006h . Each science image has an approxi- 
mate exposure time of ~0.6 h per pixel and a limiting 
magnitude of ^28 mag. 

2.3.4. Subaru Suprime-Cam 

Ground based i?-band imaging of the EGS was car- 
ried out with the Subaru Telescope as a part of the Sub- 
aru S uprime-Cam Weak-Lensing Survey (Mivazak i et al.l 
[2001 . The Suprime-Cam field of view is 0.25 deg in size, 
covering the whole IRAC map in four pointings. The 
total exposure time was 30 minutes for each pointing, 

2 http:/ /www. astro. uvic.ca/grads/gwyn/cfhtls/index. html 

3 http:/ /www4.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/cadcbin/cfht / wdbi.cgi / cfht / quick/form 

4 http://iraf.noao.edu/ 
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Fig. 1. — Layout of all the imaging surveys covering the EGS. The background grey-scale image shows the IRAC 3.6 /im mosaic. The 
black square depicts the CFHTLS pointing. The intersection of both footprints defines the region of higher photometric coverage (magnified 
on the right panel), which is referred as "main region" throughout the text (see §5$ and contains the highest resolution HST data. The 
cyan and red rectangles show the CFHT12k- BRI and SUBARU i?-band coverage, respectively. The gray and purple squares show the 
MMT-u' giz and Chandra/ ACIS pointings, respectively. The magenta circle shows the GALEX survey footprint. The green rectangles show 
the Spitzer/MlPS coverage. The orange rectangle shows VLA-20cm coverage. In the main region, magnified on the right side of the figure, 
the red tiles depict the HST/ ACS footprints (the NICMOS images overlap with this area). Slack rectangles show the MOIRCS i-S^-band 
imaging, and the green and purple lines delimit the CAHA J- and i-sT s -band surveys. 



taken in four 7.5 minute exposures, in a dithering pattern 
with ~1' spacing. The approximate limiting magnitude 
is i?=26 mag. These data were downloaded from the 
SMOKA database and reduced with the Suprime-Cam 
pipeline (SDFRED, vl.O). 

2.4. Near Infrared data 
2.4.1. HST/NICMOS 

Simultaneously to the ACS observations, the Near- 
Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) covered parallel fields in the F110W and F160W 
bands (hereafter Jno and Hiqq bands, respectively), with 
a similar exposure time and a limiting magnitude of 23.5 
and 24.2, respectively. The observations were designed to 
maximize the overlap between both HST surveys. Vir- 
tually all (58 out of the 63) NICMOS frames lie within 
the area covered by the ACS primary imaging. However, 
the smaller field of view of NICMOS (~l'xl') leads to 
a total NIR coverage of only 0.0128 deg 2 , with a 90% 
overlap with the area covered by ACS but less than a 4% 
with that covered by IRAC. 

2.4.2. Subaru MOIRCS 

In addition to the optical imaging, NIR observations 
of the EGS were also obtained with the Multi-Object In- 
fraRed Camera and Spectrograph (MOIRCS) on 5 nights 
during April - May 2006 (Pis: Fukugita, Yamada) and 
a complementary run on Jun 25, 2007. The dataset 
comprises 11 pointings covering a total of 0.09 deg 2 ori- 
ented along the original strip within 52.5° <8< 53.0° and 



completely overlapping with the HST- ACS and CFHTLS 
imaging. The median exposure time per frame is ~lh for 
an approximate limiting magnitude of K s ^23-24 mag. 
The data was reduced using dedicated scripts developed 
by the MOIRCS team involving the IRAF task MS- 
CRED, plus an addi tional de- fringing process (see e.g., 
iKajisawa etafll2009t ). 

2.4.3. Palomar and Calar Alto imaging 

Given the importance of having a continuous band 
coverag e for any kind of s tudy regarding galaxy popu- 
lations (|Wuyts et al.ll2007t ). we have incorp orated in our 
data compilation the POWIR NIR catalog (jBundv et al.l 
2006). These data were acquired between September 
2002 and October 2005 using the WIRC camera in the 
Palomar 5m telescope. The total surveyed area in the 
EGS field is 2165 arcmin 2 (0.6 deg 2 ) in the if-band and 
~l/3 of that area in the J-band. The approximate lim- 
iting magnitudes are K ~ 22.9 and J ~ 21.9. Note that 
no images were publicly available for this dataset, so we 
only use the catalogs. 

In addition, we also make use of the NIR imaging ob- 
tained by the Galaxy evolution and Young Assembly 
(GOYA 7 ) project. Two photometric campaigns were car- 
ried out to obtain NIR data o f the original Groth Strip 
o=14 /l 17 m 43 s , 5= 52°28' 41" (ICristobal-Hornillos et al.l 
120031 ) and flanking fields (jBarro et al.ll2009l ). Here we 
make use of the if-band images of the flanking fields, 
observed with the Q' instrument in the 3.5m telescope 
at Calar Alto Spanish- German Astronomical Center 
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(CAHA). The frames cover a total area of ^0.24 deg 2 
to a limiting magnitude of K s ^20.7. 

Finally, we have also included data in the J-band 
fro m the narro w band survey of Ha emitters described 
in IVillar et all (|2008f ). Three 15'xl5' pointings cen- 
tered at a=U h 17 m 31 s , <5=52°28T1", a=14H7 m 31 s , 
(5=52°28T1" and a=14H8 m 14 s , 5=52°42T5" were ob- 
served in CAHA using the H2k instrument in the 3.5m 
telescope. The combined pointings cover an area of 
0.19 deg 2 to a limiting magnitude of J — 22.9. 

2.5. Mid-to-far IR data 

2.5.1. Spitzer/IRAC 

Our sample is drawn from Spitzer near/mid-IR data 
obtained as part of the Guarante ed Time Observation s 
(GTO, PI: Fazio) and presented in lBarmbv et al.l (|2008D . 
We also included additional data from the GO program 
with ID #41023 (PI: K.Nandra). The GTO IRAC imag- 
ing data at 3.6, 4.5, 5.8, and 8.0 fj,m were obtained 
over two epochs (December 2003 and June/ July 2004). 
The dataset comprises 52 different pointings that cover 
a 2°xl0' strip with approximately the same depth. To 
achieve this homogeneous coverage and due to scheduling 
issues, the width of the mosaic is slightly variable along 
the strip, ranging from 10' to 17'. The average exposure 
time per pixel is 2.5 h (9100s) in the four channels. An 
area of 1440 arcmin 2 was observed for 190 0s, 930 arcmin 2 
for 9 100s and ~ lOOarcmin 2 for > 11500s (jBarmbv et al.l 
2008). The additional GO data are located in two strips 
of width ^3.5'flanking the original strip and covering the 
declination range 52.35°<£<53.25°. All the data were 
reduced with the general Spitzer pipeline, which pro- 
vides Basic Calibrated Data, and then mosaicked with 
Mopex using a pixel scale half of the original (~0.61 arc- 
sec pixel -1 ). The details of the image quality are dis- 
cussed in § 13.21 

2.5.2. Spitzer /MIPS 24 am and 70 am 

Complementary to the IRAC observations, MIR and 
FIR observations were also obtained with the Multi- 
band Imaging Photometer for Spitzer fMIPS: lRieke et all 
120041 ) as part of the GTO and the Far- Infrared Deep Ex- 
tragalactic Legacy Survey (FIDEL). For this paper, we 
use the whole GTO+FIDEL dataset, reduced and mo- 
saicked with the Spitzer pipeline and MOPEX+GeRT 
software. The surveyed area covers approximately the 
entire 2° x 10' strip, being slightly wider on the upper 
and lower edges. The MIPS mosaic overlaps with the 
deepest part of the IRAC observations. The mean expo- 
sure time at 24 fj,m is *~7200 s per pixel, while for the 
70 /im channel it is approximately 3800 s. The approxi- 
mate limiting fluxes are 60 /iJy and 3.5 mJy, respectively. 

2.6. Radio data 

A radio survey at 1.4 GHz (20cm) of the Northern half 
of the EGS (^50% of the IRAC mosaic) was conducted 
with the Very Large Array (VLA) in its B configuration 
during 2003-2005. The AEGIS20 survey covers 0.73 deg 2 
down to 130 /iJy beam" 1 including a smaller region of 
0.04 deg 2 with a 50 fiJy detection limit (5a). The data re- 
duction and the s ource catalog , comp rising 1123 sources, 
were presented in llvison et al.l ()2007f ). 



2.7. Keck Optical spectra 

EGS has also been the target of an exhaustive and 
unique spectroscopic follow up. As one of the DEEP2 
fields, optical multi-object spectroscopy has been carried 
out from 2003 to 2005 with the Deep Imaging M ulti- 
Object Spectrograph (DEIMOS; iFaber et all I2003D on 
the Keck II telescope. The observations cover the 
spectral range 640<A<910 nm with a resolution of 
0.14 nm. The DEEP2 DR3 6 contains spectroscopic 
redshifts for 13,867 sources at 0<z<1.4 with a median 
redshift z=0.75. The targets were selected from the 
CFHT12k BRI images (within 1.31 deg 2 ) in the mag- 
nitude range 18.5<i?<24.1. The spectroscopic redshifts 
for around 70% of the sources are labeled with a high 
quality flag (values of 3 and 4, meaning >95% success 
rate). Lower quality flags are considered unreliable and 
will be excluded from our analysis here and in Paper II. 

We complemented the DEEP2 spectroscopy of z<1.4 
galaxies with redshifts for z^ 3 s ources from the Ly - 
man Break Galaxy (LBG) survey of ISteidel et all (2003). 
This survey covers a total area of 0.38 deg 2 divided m 
several fields, one of them centered in the EGS. The ob- 
servations in EGS consist of a single 15'xl5' mask cen- 
tered at a=14 ,l 17 m 43 s , 5=52°28'48" observe d with the 
Low Resolution Imaging Spectrometer (XRIS: lOke et all 
1995) on Keck. The spectra cover the 400-700 nm 
range with a median resolution of 0.75 nm. The targets 
were pre-selected b ased on the LBG color-color criteria 
(jSteidel et al.l [1996) including only candidates brighter 
than i?=25.5. The EGS catalog contains a total of 334 
LBG candidates in the surveyed area. Out of them, 193 
are spectroscopically confirmed to be at z~3. Unfortu- 
nately, the overlap with the IRAC frame is not complete 
(and some of the galaxies are extremely faint in the IRAC 
bands), and we were only able to identify 243 (72% of the 
spectroscopic sample) LBGs in our 3.6+4.5 /xm selected 
catalog (we give more details on these sources in Paper 
II). 

3. DESCRIPTION OF THE SAMPLE SELECTION 

The dataset described in the previous Section was 
used to obtain UV-to-FIR Spectral Energy Distributions 
(SEDs) for all the sources detected in the EGS IRAC 
survey. This merged photometric catalog was built fol- 
lowing the procedure described in PG05 and PG08. Here 
we review all the basic steps of the method, emphasizing 
the improvements introduced for this paper concerning 
the extraction of the IRAC catalog and the band merging 
procedure. 

3.1. The IRAC-3.6iim+4-.5iim selection 

The source detection in the IRAC data was carried out 
separate ly in the 3.6 nm and 4.5 /im images using SEx- 
tractor (jBertin fc Arnoutsl Il9961 The complementary 
detection in the slightly shallower 4.5/xm-band helps to 
alleviate the source confusion problems arising from the 
PSF size and the remarkable depth of the IRAC data. 
Both catalogs were cross-matched using a 1" search ra- 
dius to remove repeated sources. This produces a master 
IRAC-selected catalog containing the sources detected in 
any of the two channels. Eventually, most sources are si- 
multaneously detected in both channels. 



7 http:/ /www. astro. ufl.edu/GOYA/home. html 



6 http:/ /deep. berkeley.edu/DR3/dr3. primer. html 
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TABLE 1 

Properties of the dataset in the EGS 



Band 


A cft 


m lim [AB] 


FWHM 


Area 


Surf. Dcns.(xl0 3 ) 


Source 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



Hard X-Rayf 
Soft X-Rayf 

FUV 

NUV 

u' 

u* 

B 

9 

<?' 

^606 



R.. 



I 

Z 1 



z 

Jno 

J 

#160 

•/ 

K s 

K s 

K\ 

IRAC-3.6 

IRAC-4.5 

IRAC-5.8 

IRAC-8.0 

MIPS-24 

MIPS-70 

Radio 20cmf 
R.redshiftf.... 



0.31 nm (2-10 keV) 
1.24 nm (0.5-2keV) 

153.9 nm 

231.6 nm 

362.5 nm 

381.1 nm 

439.0 nm 

481.4 nm 
486.3 nm 
591.3 nm 
625.8 nm 
651.8 nm 

660.1 nm 

769.0 nm 

781.5 nm 

813.2 nm 

807.3 nm 

887.1 nm 
907.0 nm 

1.10 fim 
1.21 fim 
1.59 fim 
1.24 /im 

2.11 /im 

2.15 fim 

2.16 fim 

3.6 fim 

4.5 fim 

5.8 fim 

8.0 fim 
23.7 A*m 
71.4 fim 

20 cm 
640-910 nm 



3.8 xl0~ 16 erg cm 2 s 
5.3 xl0 _17 erg cm 2 s 

25.6 

25.6 

26.1 

25.7 

25.7 

26.7 

26.5 

26.9 

26.3 

26.1 

25.3 

25.9 

25.3 

24.9 

26.1 

24.7 

25.3 

23.5 

22.9 

24.2 

21.9 

20.7 

23.7 

22.9 

23.9 

23.9 

22.3 

22.3 
19.5 (60 /iJy) 
15 (3.5 mjy) 
100 fijy beam 

24.1 



0.5-6" 0.67 1.102 Chandra/ ACIS; Laird et al. (20091 

0.5-4" 0.67 1.540 Chandra/ ACIS; Laird et al. (2009) 

5.5" 1.13 10.5 GALEX GTO 

5.5" 1.13 24.7 GALEX GTO 

1.0" 0.77 148.2 MMT/Megacam 

0.9" 1 152.5 CFHTLS/MegaCam 

1.2" 1.31 101.7 CFHT-12k 

1.3" 0.77 203.8 MMT/MegaCam 

0.9" 1 163.4 CFHTLS/MegaCam 

0.2" 0.197 440.7 HST/ACS 

0.8" 1 363.9 CFHTLS/MegaCam 

0.7" 1 220.0 Subaru/SuprimeCam 

1.0" 1.31 144.3 CFHT-12k 

0.8" 1 341.1 CFHTLS/MegaCam 

1.0" 0.77 275.8 MMT/MegaCam 

1.1" 1.31 117.0 CFHT-12k 

0.2" 0.197 452.1 HST/ACS 

0.8" 1 179.0 CFHTLS/Megacam 

1.2" 0.77 214.0 MMT/Megacam 

0.7" 0.0128 252.0 #ST/NICMOS 

1.0" 0.195 67.0 CAHA-f!2fc 

0.8" 0.0128 252.0 #ST/NICMOS 

1" 0.30 30.0 Palomar-WIRC; Bundv et al. (2006') 

1.5" 0.20 18.6 CAHA-fi' 

0.6" 0.09 124.0 Subaru MOIRCS 

1" 0.70 34.0 Palomar-WIRC; Bundv et al. (20061 

2.1" 0.48 315.0 Spitzer GTO 

2.1" 0.48 274.4 Spitzer GTO 

2.2" 0.48 129.7 Spitzer GTO 

2.2" 0.48 115.2 Spitzer GTO 

5" 0.79 30.0 Spitzer GTO 

19" 0.69 6.0 Spitzer GTO 

4.2" 0.73 1.538 VLA; Ivison et al. (20071 

1.31 10.343 DEEP2 



Note. — 
tData drawn from a catalog. 

Col(l) Name of the observing band. 

Col(2) Effective wavelength of the filter calculated by convolving the Vega spectrum iCol ina & Bohlin 1994|) with the transmission curve of 
the filtcr+dctcctor. 

Col(3) Limiting AB magnitude (except for the X-ray catalogs) of the image estimated as the magnitude of a SNR— 5 detection (sec ^ for 
details on the measurement of the photometric errors). 

Col(4) Median FWHM of the PSF in arcscconds measured in a large number of stars (sec 5 I5.4I I. 
Col(5): Area covered by the observations in deg 2 . 

Col(6): Source density per square degree up to the limiting magnitude given in Col(3). 
Col(7): Source from where the data were obtained. 



The average survey depth is remarkably homogeneous 
across the strip, t exp ^10ks, with the exception of two 
small areas with lower exposure at the top (5>53.525°) 
and bottom (5<52.025°) of the mosaic. We took into 
account the lower exposure times near the edges of 
the images by defining two different areas: a shal- 
lower region with exposure time shorter than 3800 sec- 
onds (N(frames)<20), and a deeper region covering the 
majority of the strip. The detection was carried out 
with different SExtractor parameters in each region, us- 
ing a more conservative configuration for the shallower 
region. Then, we used this (more restrictive) cata- 
log to purge some low significance detections in the 
other catalog within an overlapping area between them 
(N(frames)=18-25). The purged catalog restricted to the 
area with N(frames)>20 constitutes our master photo- 
metric catalog, and it covers an area of 0.50 deg 2 . 

After the detection of sources, we removed spurious 
sources in the wings of bright stars (where the PSF shows 
bright knots). For that purpose, first we made a prelim- 
inary detection of star-like sources based on the IRAC 



color-color criteria of lEisenhardt et al] ()2004 see Sec- 
tion I5.4[) . Then, we eliminated detections within the 
typical distance where the contamination from the star is 
significant (typically ~9" for the typical star magnitudes 
and depth of the EGS observations). After masking the 
regions around stars, the total area covered by the cata- 
log is 0.48 deg 2 . 

Aperture photometry for all 3.6 /xm+4.5 ^tm detected 
sources was measured with our own dedicated software 
(which takes into account pixel fractions appropriately; 
see PG08) in all four IRAC images, previously regis- 
tered to the same WCS. The flux is measured in the four 
IRAC bands simultaneously. If a source is undetected 
in the shallower bands (i.e., [5.8] and [8.0]) we still mea- 
sure an upper limit flux as 3 times the rms of the sky. 
The majority of the IRAC sources are unresolved in the 
3.6 /im image (FWHM ^2"). However, most of them 
are not point-like, but slightly extended. Consequently, 
PSF fitting is not effective and photometr y is best mea- 



sured with small circular a pertures (PG08; Barmbv et al] 
[20081 IWuvts et all [20081: lilbert et al j [2009f) . The flux 
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measurement in all bands was carried out at the posi- 
tions specified in the IRAC master catalog. We used 
a 2" radius aperture and applied aperture corrections 
for each band derived from the point spread function 
(PSF) growth curves. The values of the correction 
are [0.32±0.03,0.36±0.02,0.53±0.02,0.65±0.02] mag at 
[3.6,4.5,5.8,8.0] fj,m. The errors account for the typical 
WCS alignment uncertainties. For a small number of ex- 
tended sources (~2% of the total catalog, and 75% of 
them presenting [3.6]<22.3), the 2" aperture tend to un- 
derestimate the total magnitudes (by more than a 10%). 
These sourc es are typically bright nearby galaxies whose 
iKronl ()1980l ) radius is larger than ^4.5" . The flux mea- 
surement for these sources was performed in larger aper- 
tures enclosing the full object and applying the extended 
source aperture corrections given in the Spitzer/IRAC 
cookbook. 

The uncertainties in the IRAC photometry were com- 
puted taking into account the contributions from the sky 
emission, the readout noise, the photon counting statis- 
tics, the uncertainties in the aperture corrections, and 
a 2% uncertainty from the zero point absolute calibra- 
tion (jReach et alJl2005ft . We did not assume the uncer- 
tainties resulting from SExtractor flux measurement. In- 
stead, we used a more realistic method to determine the 
background noise that takes into account the effects of 
pixel-to-pixel signal correlation. This procedure has also 
been applied to measure the photometric uncertainties 
in all the bands and is outlined in §4.31 Nevertheless, 
a straightforward comparison to the SExtractor errors 
indicates that the noise correlation does not introduce a 
significant contribution to the flux uncertainties at bright 
magnitudes, leading to a median increment <0.02 mag 
up to [3.6,4.5]~23.75. 

A detailed comparison of our I RAC catalog o f the ECS 
survey to the one published by iBarmbv et al.l (|2008| ) is 
presented in § 13.31 including a discussion on the source 
confusion levels. 

3.2. Completeness and limiting magnitude of the IRAC 

catalog 

We estimated the completeness of the IRAC catalog 
by analyzing the recovery of simulated sources added in 
the mosaicked images. The simulations were carried out 
in the central regions of the mosaic where the coverage 
is uniform (t exp >10ks). Artificial sources spanning a 
wide range of sizes (from 1" to 6" ) and brightnesses 
([3.6]=16-25 mag) were created on the IRAC images at 
random locations. The number of simulated sources was 
chosen to be representative of the Poisson uncertainty in 
the observed number densities. The source detection and 
photometry was performed again in the simulated images 
keeping the same SExtractor parameters as in the orig- 
inal frames. The success rate recovering the simulated 
sources allow us to estimate the completeness level as a 
function of magnitude. Figure [5] summarizes the com- 
pleteness analysis in the four channels. For simplicity, 
we show only results for point-like sources. The com- 
pleteness for extended sources is typically ~10% lower 
at faint magnitudes, however, these sources represent a 
very small fraction of the catalog at these magnitudes. 
Figure [2] also shows the source density as function of the 
magnitude in each band. We find that the catalog is 85% 
complete for point sources with [3.6,4.5]=23.75 mag, and 



75% complete at [3.6,4.5]=24.75 mag. In the two other 
channels, the detection efficiency is significantly lower, 
with 85% completeness at [5.8,8.0]=22.25, and the com- 
pleteness dropping rapidly beyond that magnitude. Note 
that the forced photometric measurement in these bands 
provides a significant number of <5er (vertical dashed 
lines in Figure [5]) detections that would be missed other- 
wise. 

The lower panels of Figure [2] show the photometric un- 
certainties as a function of magnitude in the four IRAC 
bands. The red and green lines indicate the median 
and the level enclosing 90% of the distribution, respec- 
tively. For the region with the deepest coverage, we es- 
timated a 3cr limiting magnitude of [3.6,4.5]~24.75 and 
[5.8,8.0]^22.90 from the median value of the sky rms in 
our default photometric apertures (2" radius). 

The 3.6+4.5 /im catalog contains 70,048 and 99,618 
sources with [3.6]<23.75 and [3.6]<24.75 respectively. 
The median magnitude of the sample up to [3.6]<24.75 is 
[3.6]=23.14, and 75% of the sources present [3.6]<23.96. 
Note that these numbers correspond to the IRAC 
3.6+4.5 /zm catalog before applying the de-blending tech- 
nique discussed in § [5] Consequently, the number of 
sources quoted above is lower than in the final catalog 
(see Table [2]). 

3.3. Comparison to Barmby et al. IRAC selected catalog 

Here we compare o ur IRAC pho t ometr ic catalog 
to that published by IBarmbv et al.1 ({20081 hereafter 
BAR08). BAR08 used the same dataset (except for the 
GO flanking regions) and obtained final mosaics in all 4 
IRAC bands with a very similar reduction to ours. Con- 
cerning the source detection, our method is slightly dif- 
ferent from BAR08 since we detect galaxies in both the 
3.6 and 4.5 fim channel instead of only using the bluer 
band. Our dual detection technique helps to alleviate the 
source confusion problems arising from the PSF size and 
the remarkable depth of the IRAC data. In addition, we 
measure fluxes in the four channels simultaneously, ob- 
taining upper limit values for undetected sources in the 
shallower bands ([5.8] and [8.0]). Note also that we have 
increased the resolution of our catalog by deconvolving 
blended sources using higher resolution information from 
ground-based observations f ij4.ll) . However, for the sake 
of clarity, we compare here the BAR08 catalog with ours 
before carrying out the deblending procedure. 

BAR08 measured aperture photometry with SEx- 
tractor. The publicly available catalog 1 includes 
MAG_AUTO and MAGJSO measurements, jointly with 
aperture magnitudes for several radii corrected to total 
magnitudes with empirical PSF corrections. We compare 
our photometry to the magnitudes measured in the 3.5 
pixel aperture (2.1") by BAR08. Their aperture correc- 
tions agree with our measurements for the 4" diameter 
aperture within the errors (due to alignment uncertain- 
ties). 

BAR08 and our catalog are cross-correlated in the re- 
gion of highest exposure (t e3; p>4ks) using a 1" radius. 
The comparison of the WCS between the two mosaics is 
in very good agreement, with a rms of <0.05" . We find 
41,514 and 52,130 sources in common up to [3.6]<23.75 
and [3.6] < 24. 75, respectively. 

1 http:/ /www. cfa.harvard.edu/irac/egs/ 
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[5.8] [8.0] 



Fig. 2. — The upper panel of each quadrant shows the histogram of IRAC magnitudes for the sources in our sample selected in the 
IRAC-3.6/xm+4.5/im data. The histograms are shown up to the 2a limiting magnitude in each band. The blacks continuous curve depicts 
the detection efficiency as estimated from recovery rate of simulated punctual sources. The vertical dashed lines represent 5 times the 
median sky-rms (5<r) measured in a large number of 2" radius apertures. The lower panels in each quadrant show the distribution of 
photometric uncertainties as a function of the magnitude in each IRAC band. The red and green lines indicate the median and 90% of the 
error distribution as a function of magnitude, respectively. The horizontal dashed line shows the SNR~5 limit; the intersection with the 
red line indicates the value quoted in Table [l] for the IRAC bands. Note that these values are slightly lower than the 5cr sky-rms. 
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Attending to the density of sources per unit area, we 
find that our catalog includes 11%±5% and 18%±7% 
more sources than the catalog published by BAR08 at 
[3.6]<23.75 mag and [3.6]<24.75 mag, respectively. The 
uncertainties in these measurements were estimated by 
comparing number counts in 0.5 mag bins, including 
Monte-Carlo simulations on the photometric errors. 

The different source densities are a consequence of their 
more conservative SExtractor detection threshold. We 
have carefully chosen the SExtractor parameters differ- 
entiating between the regions with high/low coverage 
( H3.1[) . trying to push down the detection limits as much 
as possible without degrading the reliability of the en- 
tire catalog. As a consequence, our catalog recovers a 
larger number of sources at fainter magnitudes, and the 
completeness of our catalog is larger than BAR08 for the 
same magnitude. For example, they quoted a < 50% 
completeness for point-like sources at [3.6]=23.75 mag, 
compared to our estimated 85%. Unfortunately, lower- 
ing the detection threshold inevitably increases the num- 
ber of spurious detections. However, in the context of a 
merged multi-band photometric catalog, spurious detec- 
tions can be efficiently identified as sources only detected 
by 1RAC (cf. N(band)<5), and we will show that the re- 
liability is > 97%, with false detections located almost 
uniquely around very bright sources (see § 15. 3[) . 

Figure[3]shows the comparison of the photometric mag- 
nitudes in the four bands for both catalogs. In both 
cases the flux was measured on circular apertures and 
corrected to total magnitudes using aperture corrections. 
We have corrected the comparison by a constant value 
of -0.05 mag in [3.6] and [4.5], and by -0.04 mag and 
-0.03 mag in [5.8] and [8.0], respectively. Such small off- 
sets are attributed to slight differences in the data re- 
duction (final absolute calibration, frame stacking, regis- 
tering and mosaicking) and in the aperture corrections. 
Despite the small offsets, the overall results in the four 
IRAC bands are in good agreement. Note that the aver- 
age photometric error in our catalog for each magnitude 
bin (green bars) encloses lex of the values around the 
median difference (red and cyan lines). 

The lower panel of each plot in Figure [3] shows the 
comparison of the photometric errors in BAR08 and in 
our catalog. Our quoted photometric uncertainties tend 
to differ from BAR08, specially at faint magnitudes. In 
contrast with that paper, we have considered zero-point 
and WCS uncertainties, resulting on slightly larger un- 
certainties (^0.05 mag) in our catalog for bright sources 
up to [3.6][4.5]~21 mag and [5.8][8.0]~19 mag. At fainter 
magnitudes, the photometric uncertainties increase with 
magnitude at a faster rate in the catalog of BAR08. The 
cause for this difference is the procedure to measure the 
background noise. Similarly to BAR08, we estimated 
this value from the sky variance measured in circular 
apertures at different locations of the images that are 
empty of sources ( ^4.3p . However, the definition of an 
empty region depends on the limits of source detection. 
Therefore, given our higher detection fraction, our sky 
regions would contain, in principle, lower signal pixels 
effectively decreasing the rms. In addition, our photo- 
metric procedure estimates uncertainties on a source-by- 
source basis studying the background around each object 
in a independent way, while BAR08 relied on SExtrac- 
tor photometric errors and applied a correction to them 



based on the average properties of the mosaic. Neverthe- 
less, our estimates of the photometric errors are consis- 
tent with the observed scatter of the comparison between 
our photometry and that measured by BAR08. 

4. MULTI- WAVELENGTH PHOTOMETRY: THE RAINBOW 
CATALOG 

Using the whole dataset available in the EGS field, we 
created a multi-wavelength photometric catalog for the 
IRAC selected sample described in the previous Section. 
For that purpose, we used the Rainbow software package, 
described in detail in PG05 and PG08. This software was 
created to: (1) cross-correlate multi-band catalogs and 
obtain consistent (aperture matched) photometry on the 
different bands to build a UV-to-FIR SED; (2) estimate 
stellar parameters, such as photometric redshifts, stel- 
lar masses, and SFRs from those SEDs. In the rest of 
this Section we will describe the photometry procedure, 
and in Paper II we will present the methods to estimate 
photometric redshifts, stellar masses, and SFRs out the 
SEDs. 

4.1. Cross-correlation and source deblending 

The Rainbow code starts from a primary selection cat- 
alog (in our case, IRAC selected) and obtains merged 
photometry and spectroscopy in other bands. The first 
step is to identify the counterparts of the IRAC sources 
in the other bands, where SExtractor catalogs have been 
built following typical procedures. These catalogs are 
cross-correlated to the 3.6+4.5 /xm positions using a 2" 
search radius. An exception to this rule are the MIPS, 
radio and X-ray catalogs. For the MIPS and radio bands, 
we used a 2.5" and 3" radius, respectively, recognizing 
possible alignment and center estimation problems of the 
order of one pixel. For the latter, instead of using the 
WCS of the X-ray sources, we cross-matche d to the po- 
sition s of the IRAC counterparts (given in lLaird et al.l 
l2009h using a 1" radius. These authors used as reference 
the IRAC catalog of BAR08, which cover a slightly lower 
area than ours. Thus, for sources outside of the BAR08 
mosaic we cross-matched to the X-ray coordinates using 
a 2" radius, (see § 15.21 for more details). 

The IRAC sources are identified in this way with ob- 
jects in all the other catalogs. One of the optical catalogs 
(typically the deepest; in the EGS, the Subaru f?-band 
data) is used as a reference to narrow the following cross- 
correlations and alleviate confusion problems present in 
the IRAC images. As the cross-correlation to the cata- 
log of spectroscopic redshifts is done to the coordinates 
of the counterpart in the reference (optical) band, it is 
possible to choose a more reliable 0.75" search radius for 
this catalog. 

Before the cross-matching procedure is carried out 
between a pair of images, these are re- aligned locally 
within a 4' x4' square region using the positions of several 
sources (typically more than 20) as reference. The mean 
rms between the central coordinates of matched sources 
in optical/NIR images is typically <0.1" and <0.2" be- 
tween the IRAC and the ground-based images. This pro- 
cedure allows to overcome small misalignment problems 
between the frames, and assures a reliable identification 
of counterparts, and an accurate positioning of the pho- 
tometric aperture in all bands. It also allows to obtain 
reliable photometry (in the appropriate aperture) even if 
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[5.8] [8.0] 

Fig. 3. — Comparison of the observed magnitudes (upper panels) and the photometric errors (lower panels) in the four IRAC bands for 
the sources in common between our IRAC-3.6+4.5 ^m catalog and BAR08 catalog. The cyan line shows the median value of the magnitude 
difference as a function of magnitude. We have corrected the comparisons by a constant value of -0.05 mag in [3.6] and [4.5], and by 
-0.04 mag and -0.03 mag in [5.8] and [8.0], respectively. Such small offsets can be attributed to slight differences in the data reduction and 
aperture corrections. The green bars indicate the average photometric errors per magnitude bin in our catalog. The red lines enclose la of 
the distribution centered in the median value (cyan line) . The vertical dashed line indicate the 5cr limiting magnitude in our catalog. The 
photometry in the four bands is consistent up to the ~ 85% completeness limit. The uncertainties in BAR08 arc 5% to 10% larger at faint 
magnitudes, probably as a result of the slightly different procedure applied to measure the sky background. 
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Fig. 4. — Example of the deconvolution procedure for unre- 
solved sources in IRAC images. After the deblending, the unique 
irac source (iracl64074) becomes two separate sources (iracl64074 
.1 and .2). Top-left (A): IRAC-3.6^m image (0.61"/pixel) showing 
a blended source. The red circle depicts a 2" radius aperture (our 
default aperture for isolated sources). The green circles show 0.9" 
apertures centered at the positions of sources in the reference im- 
age. Top-right (B): HST/ ACS reference image (0.02" /px) after 
applying a 4' x4' local WCS re-alignment, showing the individual 
sources. Bottom-left ( C): Model of the blended source obtained by 
convolving the PSF of HST /ACS to the PSF of IRAC-3.6/mi in 
(B), registering it to image (A) and scaling each source to the flux 
in the 0.9" apertures (green) in (A). Bottom-right (D) Residual 
from the subtraction of the model and science images (C-A). 

the source is very faint and/or undetected in an individ- 
ual image. 

The combination of the remarkable depth and the ^2" 
FHWM of the IRAC observations inevitably leads to is- 
sues of source confusion, specially around crowded envi- 
ronments. However, based on the (ground-based) ref- 
erence image, it is possible to deblend IRAC sources 
which have not been separated by SExtractor in the 
original IRAC images and lie at least 1" away (half 
the FWHM, chosen as our resolution criterion). HST 
images reveal that the multiplicity is larger, but the 
deconvolution of sources separated by less than 1" is 
very uncertain. When multiple counterparts are found 
in the optical/NIR images during the cross-matching, 
the IRAC photometry is recomputed foll owing a decon- 
voluti o n method sim il ar to t hat used in iGrazian et al.l 
<l200?lh. iWuvts et all (pOOSh . iWilliams et all lj2009h or 
iWang et all (|2010l ). 

In our case, first, the coordinates of the photometric 
aperture are re-positioned to that of the optical/NIR 
counterparts. Then, the PSF of the higher resolution 
image is convolved to the IRAC PSF, and the flux of 
each source is scaled to match that of the real IRAC 
sources measured in 0.9" apertures (after re-centering to 
the positions of the optical counterparts). Finally, total 
magnitudes are computed applying an aperture correc- 
tion of [1.30±0.07,1.02±0.08,1.2±0.10,1.44±0.14] mag in 
the [3.6,4.5,5.8,8.0] fjjm bands, respectively. Figure 0] il- 
lustrates the deconvolution procedure using a HST/ ACS 



image as reference. The red and green apertures depict 
the standard 2" aperture (for isolated sources) and the 
0.9" aperture, respectively. Pixel-by-pixel variations in 
the residual from subtracting the model PSF do not ex- 
ceed a 5% within the 0.9" and 2" apertures. We also 
checked that the average rms (~3%) is well within the 
photometric error of the sources. The analysis of an aver- 
age PSF, derived from observed sources across the image, 
indicates that for the typical separation between blended 
sources, ~2.2" (>1.8" for 75% of them), the flux contam- 
ination from the nearby neighbor does not exceed a 10% 
for sources with a flux ratio around 1:2-3. Approximately 
75% of the blended sources present flux ratios lower than 
1:3.5. 

After applying the deblending method, our IRAC 
selected catalog contains 76,936 (113,023) sources to 
[3.6]<23.75 (24.75). This means that we were able to de- 
blend 8% of the sources in the original IRAC catalog built 
with Sextractor (presented in § 13. 1|) . and 16% of the final 
catalog of 76,936 sources were deblended (typically, each 
blended sources was a combination of 2 sources) We find 
no significant difference in the brightness distribution of 
the blended sources compared to the resolved sources. 

In the following sections, and in Paper II, we 
will analyze the SEDs and physical properties of 
the IRAC sources, concentrating in the sample with 
[3.6]<23.75 mag, which count with more accurate IRAC 
photometry (SNR>8). Therefore, this will be the work- 
ing sample for the rest of the paper unless explicitly 
stated otherwise. Nonetheless, all the procedures dis- 
cussed in the following are also applied to the sources 
down to [3.6]<24.75 mag. Although these objects are 
not included in the accompanying catalog (presented in 
§ El restricted to the most reliable detections) , it is pos- 
sible to retrieve their data through our online database 
(see SO]) . 

4.2. Merged photometry 

The photometry is carried out in the same (Kron) el- 
liptical aperture in all bands. The parameters of that 
aperture are obtained from a reference image (the same 
as for the cross-matching procedure) whose resolution is 
representative of the entire dataset. Normally, this refer- 
ence image is a ground-based optical/NIR frame with a 
PSF of approximately 1-1.5", which is easy to translate 
to other ground-based images avoiding aperture issues. 

The bands are sorted according to depth to facilitate 
the cross-match to the optical bands. The typical aper- 
ture band for the majority of the sources (83% of the 
sample) is the Subaru i?-band, followed by the CFHTLS 
«'-band (4%), MMT -i-z (3% each) and MOIRC-X 
(2%). For the remaining 5% of the sources, the aper- 
ture is computed from other bands (CFHTLS-r, isi4, 
V505) that account for less than 1% of the total each. 
The Subaru imaging was preferred to the CFHTLS as 
primary aperture band due to the uniform coverage of 
the whole IRAC mosaic (the CFHTLS frames cover only 
~60% of the IRAC survey). In order to ensure that the 
aperture radius is always large enough to enclose the full 
PSF profile in all ground-based images, we established 
a minimum aperture radius equal to the worse value of 
the seeing (1.5"). Once the best photometric aperture is 
defined, if no counterpart is found in one band but data 
exist at that position, the flux is measured within the 




Fig. 5. — Multi-band images of an IRAC selected source 
(irac096525) in the [3.6], R, I and ign bands (from left to right and 
top to bottom). The size of the frames is 15" xl5" . In the three 
optical bands, the flux is measured in elliptical apertures (red) 
whose parameters are determined in the R-band image. Note that, 
despite the different image resolutions (and multiple HST counter- 
parts), the accurate WCS re-centering (rms<0.15") allows to place 
the aperture correctly, recovering the flux of all sources. For the 
IRAC bands, we use fixed circular apertures of 2" radius (blue) 
and we apply a correction to the total magnitude. 

same aperture. Note that the local WCS re-centering of 
the images allows an accurate positioning of that aper- 
ture. In this way we recover fluxes for very faint sources 
whose detection was missed by SExtractor. When the 
counter-part source is too faint to be detected, the sky- 
rms value is stored to be used as an upper limit in the 
SED analysis. 

Some bands are deliberately excluded from the merged 
photometric measurement described above due to a sig- 
nificant difference of the resolution compared with the 
reference band or because images are not available. For 
all these bands, the photometry is obtained differently 
and later incorporated in the merged catalog during the 
cross-correlation procedure. For the four IRAC chan- 
nels, the optical/NIR Kron apertures are typically not 
large enough to enclose the entire PSF. Therefore, we 
keep the values measured in §3.11 with aperture photom- 
etry. Given the comparatively large PSFs of the MIPS 
(24 (im, 70 (Jim) and GALEX (FUV, NUV) images, the 
flux measurement in these bands was carried with IRAF- 
DAOPHOT and SExtractor, respectively (see PG05 and 
PG08). Figure [S] illustrates the different resolutions and 
aperture sizes involved in the photometric measurement 
(see also the captions in Figure l9|10|lll and H"2"j) . 

The HST- ACS bands were included in the general pho- 
tometric procedure. Although the much higher resolu- 
tion may lead to multiple cross identifications even when 
compared with the optical ground based images, the high 
spatial resolution (after the local re-alignment described 
above) allows a reliable cross-matching within 0.15" (rms 
of the local WCS solution for HST images). As shown 
in Figure [5l apertures were placed correctly even when 
multiple counterparts are identified, and the photometry 



includes the fluxes from all them. We have conducted an 
additional test to check the accuracy of the ACS photom- 
etry measured with this method. We compared the ob- 
served photometry to synthetic magnitudes derived from 
SED templates fitted to the multi-band photometry (ex- 
cept ACS) of spectroscopically confirmed galaxies. We 
find a very small offset (< 0.02 mag) and a scatter con- 
sistent with the typical photometric errors in the bands 
(see §3.1.3 of Paper II for more details). 

T he NIR source cata log of the Palomar/WIRC sur- 
vey (|Bundv et al.l [20061 is also included in the merged 
photometric catalog, although no images are available to 
match apertures. In this case, we use the SExtractor 
MAG_AUTO value of the closest neighbor in the J- and 
K- bands. 

4.3. Photometric uncertainties 

The uncertainties in the photometry are derived si- 
multaneously to the process of flux measurement in each 
individual band. As mentioned in H3.ll the photometric 
errors obtained using SExtractor often underestimate the 
true backgr ound noise due to signal correlation on adja- 
cent pixels (|Labbe et al.ll2003l iGawiser et al J 12006( 1. In 
order to properly account for this effect, we estimated the 
flux uncertainties in three different ways, as described in 
PG08 (Appendix A). First, we measured the background 
noise in a circular corona of 5 " width around each source. 
This procedure is similar to that used by SExtractor and 

1/2 

provides uncertainties ctap oc N p ' ixels . In addition, we 
measured the average sky signal on non-connected ar- 
tificial apertures built with random pixels around each 
source, with the same size as the photometric aperture, 
and containing only pure sky pixels (rejecting pixels >5cr 
the rms obtained with the first method). Finally, we also 
estimated the back ground noise following the method by 
lLabbe et al.l ((2003) ■ The flux measured on several pho- 
tometric apertures around the source, identical to the 
one employed for the photometry, is fitted to a Gaussian 
function to yield the rms background fluctuation. The 
sky background is set to the resulting average value of 
the three methods, and the final photometric uncertainty 
is set to the largest estimate, typically, the one measured 
with the second method. 

5. PROPERTIES OF THE MERGED PHOTOMETRIC 
CATALOG 

The exposure time of the IRAC mosaics in the four 
channels is remarkably homogeneous across the mosaic. 
However, the coverage of the strip at other wavelengths is 
patchy and discontinuous (see Figure [T]). For the remain- 
ing of this Section and in Paper II, we will differentiate 
between the region of the IRAC image covered by the 
CFHTLS, and the rest. The CFHTLS/IRAC common 
region (52.16° < 5 <53.20° & 214.04°< a <215.74°) 
has been also surveyed with Subaru, the CFHT12k in- 
strument, MMT/Megacam, HST/ACS, ffST/NICMOS, 
and GALEX , so it constitutes the zone with the highest 
data quality (hereafter, main region), where the SEDs 
are sampled with the highest band coverage. The area of 
this region is 0.35 deg 2 (~68% of the total) and contains 
53,030 (76,936) sources down to [3.6]<23.75 (24.75). The 
0.13 deg 2 outside the main region (hereafter, flanking re- 
gions) , also present a solid SED coverage, relying mostly 
in the MMT and SUBARU bands. Nevertheless, it 
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Fig. 6. — The central panel in each plot shows the color magnitude diagrams in the u* , R, i' and K bands with respect to [3.6] (top to 
bottom, left to right). The density map in grey scale shows the distribution of IRAC sources detected in each band. For sources missing 
(not detected by Sextractor) in a given band but detected in any other, we force a flux measurement at the position of the existing source 
using the same Kron aperture. The green dots depict sources for which we are able to recover a positive flux in this forced measurement. 
The blue dots depict the value of the sky rms in the apertures where the forced measurement failed (i.e., the integrated flux was negative). 
The cyan triangles are sources undetected in any other band but IRAC. For these sources, we set an upper limiting magnitude ~0.5-l mag 
fainter than the typical magnitude for a source with SNR=2 (black line). The magenta solid line indicates the median magnitude of the 
sources with SNR~5. The red dots depict galaxies with spectroscopic redshift z>2.5. The red ellipses show the median and quartiles of 
the color-magnitude distribution for galaxies in different bins of redshift within 0<z<1.5. The upper panel of each quadrant shows the [3.6] 
brightness distribution of undetected sources (blue and cyan markers in the central panel) in each band. The right panel of each quadrant 
shows the color distribution of detected sources (histogram, median, and quartiles) in three bins of the [3.6] magnitude. 
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TABLE 2 

Number of sources in the IRAC-3.6+4.5/iM 

CATALOG 





[3.6]<23.75 a 


[3.6]<24.75 6 


Prior to deblcnding 


70,048 


99,618 


After deblending 


76,936^ 


113,023 t2 


Main c / Flanking regions 


53,030/23,906 


77,607/35,416 



Note. — 

a, 85% completeness level. 

b, 3<r limiting magnitude. 

C, The main region is defined as: 52.16° < 5 <53.20° & 214.04° < 
a <215.74°. The Flanking regions consist of the remaining area. 
1 1 , Only these sources have been included in the final catalog pre- 
sented in §[6] 

t2, A larger, although less complete, sample including all the sources 
down to [3.6]< 24.75 can be accessed through the web interface Ram- 
bow Navigator f^ fo~3l . 

lacks some of the best quality data (taken with ACS 
and MOIRCS), and the fraction of surveyed area by 
the WIRC- J K bands is -15% lower. There are 23,906 
(35,416) sources in the flanking regions with [3.6]<23.75 
(24.75). Table [5] summarizes the number of sources in 
the catalog according to different limiting magnitudes 
and area constraints. 

Note that the bulk of the optical SED coverage is based 
on CFHTLS and MMT data (the CFHT12k-BRI data 
are —2 mag shallower). On average numbers, their re- 
spective filter sets cover a similar wavelength range with 
comparable data quality (see Tables [1] and @|. Also, the 
subtle differences in the shape of the u and z band filters 
of each survey improve the quality of the SED coverage 
for the sources in the main region. The most remark- 
able difference between the two surveys is found in the 
homogeneity of data; whereas the CFHTLS present a 
nearly uniform data quality across the mosaic, the MMT 
data show larger variations between the four pointings 
(e.g., Ag(5o-)=26 -25 mag, AFWHM= 1.0 - 1.6" ; see 
iZhao et all [20091 for more details). Unfortunately, the 
lowest quality pointings (shallower data and highest see- 
ing) are precisely those covering the flanking regions 
(mainly the north region, S >53.20°). Furthermore, a 
small area in the north and south regions of the MMT 
mosaic (<10% of the total) is not observed in the g and 
u bands, respectively. Thus, comparatively, the average 
data quality in the flanking regions is slightly lower than 
in the main region. 

5.1. Multi-band detection and color properties. 

Table [4] shows the fraction of the entire IRAC catalog 
detected in the different bands compiled for this work. 
We give specific values for different sub-samples divided 
as a function of brightness in the selection band, detec- 
tion level (SNR), and location in the IRAC mosaic. The 
majority of sources are detected in the 4 IRAC chan- 
nels. However, while most sources (>90%) have reliable 
detections at 3.6 and 4.5 /im, only half of the sample 
is detected with SNR>5 at the two longer wavelengths. 
Note that the photometry at 5.8 and 8.0 /zm was not 
performed independently (see £13.11) . 

The fraction of IRAC sources detected in the op- 
tical bands is also elevated (—95%) for the majority 
of the bands. However, cutting the measurements at 
SNR>5, the efficiency decreases to —85% for the deep- 
est CFHTLS, MMT,' Subaru and ACS bands, and to 
50%-60% for the u and z bands (in both CFHTLS and 
MMT data) and the shallower CFHT BRI bands. For 



the GALEX data, the detection is typically lower than 
10% and 25% in the FUV and NUV bands, respectively. 
Around 12% and 6% of the sources are matched to high 
quality spectroscopic redshift estimates in the main and 
flanking regions, respectively. The lower efficiency in the 
flanking regions is caused by the inhomogeneous DEEP2 
coverage of the IRAC mosaic. 

The NIR coverage of the strip surveyed by IRAC is 
also discontinuous. The WIRC catalog provides the most 
uniform coverage, including 40% and —100% of the area 
in the main region in the J and K bands, respectively. 
The fraction of IRAC sources detected in each band is 
very similar, -50-60% (-20% with SNR>5). The deep- 
est NIR observations are those taken in CAHA-J and 
Subaru-MOIRCS-if s bands. The latter covers -1/4 of 
the main region up to 1 mag deeper than the WIRC-isT 
data, presenting a much higher detection fraction, —90% 
(70% with SNR>5). The CAHA-J data covers 40% of 
the main region recovering —30% (SNR>5) of the IRAC 
detections, a slightly larger fraction than WIRC- J. 

We note that for most bands, the fraction of sources 
detected with SNR<5 can be significantly higher than 
the overall value. In the shallowest bands, this is caused 
by the higher photometric uncertainties (e.g., the WIRC 
J and K bands). However, for most bands this is the re- 
sult of the enhanced detection achieved with the forced 
photometric measurement. For each IRAC source unde- 
tected in a given image, but having a counterpart in any 
other band, we still measure the flux in the same aper- 
ture at the position of the counterpart. With this method 
we increase the detection efficiency at fainter magnitudes 
recovering sources that would be missed otherwise. How- 
ever, given the extreme faintness of these sources, some 
of the measurements produce low significance detections 
(SNR<2). As a precaution, to preserve the overall qual- 
ity of the SED, we do not include these values (photomet- 
ric uncertainties >0.4 mag) in the SED fitting procedure, 
nor in the estimate of the physical parameters in Paper 
II. Nevertheless, we keep these values as they can be 
useful as upper limits. Another estimate of the limiting 
magnitude in each band is obtained from the the value 
of the sky rms in failed forced measurements (i.e., those 
for which the measured flux in the aperture is negative) . 

Figure [5] illustrates the characteristics of the different 
flux measurements in four bands probing representative 
wavelengths in the UV-to-NIR range. The central panel 
of each plot shows a color-magnitude diagram in [3.6] vs 
u*Ri' and K, respectively. The grey-scale density map 
depicts the distribution of standard (non-forced) photo- 
metric measurements, typically detected with SNR>5 
(magenta line). Forced detections for which we obtain a 
valid flux or a sky rms value are shown as green and blue 
dots, respectively. The cyan markers are sources detected 
in IRAC only (for which we assign arbitrary magnitudes 
in each band). For the deep R and i' bands, which 
present the highest detection efficiency, forced measure- 
ments account for less than 5% of the flux measurements, 
and present a median SNR<2 (black line). In the u* 
band, the source density for a similar limiting magnitude 
is much lower, and consequently the fraction of forced de- 
tections is higher, around 20%. 

The upper panels in each quadrant of Figure [5] shows 
the fraction of undetected sources in each band. In the R 
and i' bands these sources make up for <4% of the total 
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Fig. 7. — The central panel in each plot depicts the source density as function of the number of bands in which the source was detected 
(a positive flux was measured), N(band). and the magnitude in the [3.6] channel, for the main region (left) and the flanking regions (right). 
The color contours contain (from the inside out) approximately 25%, 50%, 75% and 90% of the sources. The green dots mark galaxies with 
a spectroscopic redshift. The black lines indicate the 85% completeness limit of the catalog ([3.6] <23.75, vertical) and the band coverage 
limit N(band)=5 (horizontal), respectively. The histograms on top show the brightness distribution for the full sample (black) and for 
sources with N(band)<5 (blue). The histograms to the right show the N(band) distribution of the full sample (black line), IRAC faint 
sources ([3.6]>23.75; red area) and galaxies with spectroscopic redshifts (green area). 



sample. Most of them are IRAC-only sources (plotted in 
cyan), whereas the rest are at least marginally detected 
in one other (typically red) band, but the forced mea- 
surement fails for the particular band shown in the plot 
(blue). The fraction of undetected sources in the u* band 
(and other shallower optical bands, e.g., CHFT12K-BRI) 
increases to 10%, being in this case dominated by failed 
forced measurements at the positions of R or i' detec- 
tions. These kind of measurement are also predominant 
for undetected sources in the if -band. 

Note that some of the IRAC-only sources are relatively 
bright, 21<[3.6]<22. Most of them are detected in deep 
if -band observations. However, there are a few galaxies 
detected only in the IRAC bands (typically the faintest 
ones). The nature of these interesting sources, poten- 
tial candidates to massive high redshift galaxies, will be 
explored in a future paper. 

Figure [6] also shows the color distribution of IRAC 
sources in bins of [3.6] mag (right panels in each quad- 
rant), and spectroscopically confirmed galaxies in sev- 
eral redshift bins (red ellipses and dots in the cen- 
tral panels). The general trend with color in the R 
and i' bands is that faint IRAC sources are on aver- 
age bluer than brighter ones. The median colors for 
the faintest [3.6] bin, 22.50<[3.6]<23.75, are R-[3.6]=1.6 
and %'— [3.6]=1.4. These colors are similar to those of 
galaxies at l<z<1.5. This is consistent with the fact 
that the median redshift for the IRAC-3.6+4.5 /iin mag- 
nitude limited sample ([3.6]<23.75) is z~l (See Paper 
II and PG08). Note also that most of the forced detec- 
tions in the R-band would qualify as IRAC extremely red 
objects (by the criteria of lYan et al]|2004 i?-[3.6]>4.0), 
that targets dusty starbursts or passively evolving galax- 
ies at z>1.5. 

In general, all galaxies tend to become redder in 
optical— IRAC colors with increasing redshifts. Indeed, 
for a typical galaxy SED, the observed optical bands 
shift into the (fainter) UV whereas the [3.6] mag becomes 
brighter as it approaches the stellar bump at 1.6 [im rest- 



frame (for z<1.2). Therefore, it is not surprising that the 
u* — [3.6] median color is redder [u* — [3.6]=2.1, typical of 
a z~l galaxy) than the median color involving the Ri'K 
bands. Note also that the observed UV progressively 
shifts into the Lyman break, producing the large frac- 
tion of z>2.5 (red dots in the top-left quadrant of Fig- 
ure galaxies that are u-dropouts. Finally, we find that 
the if —[3.6] color presents an almost constant value as a 
function of [3.6] mag for a given redshift range, becom- 
ing progressively redder as we move to higher redshifts: 
if-[3.6]<0 at z<0.5 and if-[3.6]>0 at z>l. Again, this 
is consistent with the fact that both if and [3.6] transit 
through the peak of the stellar bump for z<l, cha nging 
their relative positions (see, e.g.. lHuang et al.ll2004f) . An 
interesting consequence of the red if —[3.6] colors of high- 
z galaxies is that IRAC-3.6 /im observations are equiva- 
lent, in terms of source densities, to a if-selected sample 
down to slightly deeper limiting magnitudes. 

Figure [7] shows the density map (central panel) and 
histograms (right panel) of the number of photometric 
bands with measured fluxes, N(band), in the main (left) 
and flanking regions (right). Typically, the SED of a 
galaxy in the main region has a median of 19 photometric 
data points. The average spectral coverage is ~8 bands 
larger than in the flanking regions, mainly due to the 
lack of data from the CFHTLS and NIR surveys. Galax- 
ies with an available spectroscopic redshift (i?<24 mag; 
green dots) typically present the highest band coverage 
(~22 and ~14 bands in the main and flanking regions, re- 
spectively), since they are relatively bright in the optical. 
Interestingly, the faintest IRAC sources, [3.6]>23.75 mag 
(red histogram), present a relatively high band coverage, 
N(band)=17 in the main region. In most cases, these are 
low significance detections in the deepest optical/NIR 
bands where the forced photometric measurement (see 
§ 14. 2|) is able to recover a flux (see, e.g., Figured]). On 
the opposite side, there is a non-negligible number of 
IRAC-faint but optically-bright sources with more than 
19 photometric data points. These constitute a popula- 



16 



TABLE 3 

Detection efficiency of X-ray, FIR and Radio sources 



Source 


Area 




Fraction 


Total 


(1) 


(2) 




(3) 


(4) 


Chandra/ ACIS 1 " 1 


0.67 de^ 


? 


0.70 


1023 (77%) 


MIPS-24 /jm (f>60/iJy)t2 


0.53 dot 




0.59 


10771 


MIPS- 70 nm (f>3500AtJy)t2 


0.50 de f 


;2 


0.61 


868 


VLA 20 cmta 


0.73 de f 


;2 


0.46 


590 (52%) 



Note. — 

fi Sources drawn from Laird et al. I 2009). These include, 815 high reliability 
identifications (based on a previous match to the IRAC catalog of BAR08) 
and 176 previously un-idontificd sources. 

f2 The area of the survey refer to the GTO + FIDBL observations that re- 
duced for this paper. 

f 3 Sources drawn from llvisoii et aU HUUTt) . 

(1) Name of the band. (2) Total area of the survey. (3) Fraction of the survey 
area overlapping with the IRAC mosaic. (4) Number of sources with IRAC 
counterparts ( [3.6] < 23. 75) and fraction of recovered sources from the whole 
catalog. 

tion of blue dwarf galaxies at intermediate redshift easily 
detected in the optical but with faint IRAC counterparts 
(i.e., not very massive). 

The upper panel of Figure [7] shows the [3.6] mag- 
nitude distribution for the full sample (black) and for 
sources with poor spectral coverage (N(band)<5; blue 
histogram). The latter defines a clearly isolated group 
in the density contours (central panel) and in the his- 
tograms (right panel) of Figure [7] These sources repre- 
sent less than 3% (4%) of the sample up to [3.6]<23.75 
(24.75). Most of them are clear IRAC-only detections 
(as the ones discussed in Figure |6]) relatively bright in 
[3.6,4.5]. However, some of the faintest sources can be 
affected by contamination of spurious sources (we give 
more details on these sources in § 15. 3[) . 

5.2. FIR, X-Ray and radio counterparts 

Tabic [4] shows the fraction of IRAC sources in our sam- 
ple ([3.6]<23.75) detected in the X-ray, FIR and radio 
surveys within the overlapping area. About 30% (20% 
with SNR>5 detections) of the IRAC sources are de- 
tected at 24 /mi (SNR=5 reached at -60 /iJy), 10% (2% 
with SNR>5) at 70 /mi (SNR=5 reached at -3.5 mJy), 
1% in the X-rays and <1% at 20 cm. For the radio 
and X-ray surveys, we c r oss-m atche d our sampl e to th e 
catalogs of llvison et "aU (|2007t ) and lLaird et~al] (|2009ft . 
respectively. As explained in § 21 the cross-match to the 
radio catalog was performed using a 3" radius, whereas 
for the X-ray catalog we followed a two step identifica- 
tion. First, we divided the X-ray sample attending to the 
availability of pre-identified IRAC counterparts (Laird 
et al. used the IRAC catalog of BAR08 as reference). 
For the galaxies identified in IRAC, we used a 1" cross- 
matching radius, whereas for the other we used a more 
conservative 2" radius. With this procedure we found 
848 out of the 882 galaxies with IRAC counterparts in 
their ca talog (815 among the 830 with high reliability 
flag; see lLaird et alJ[200a for more details). The remain- 
ing galaxies were either outside of the N(frame)>20 re- 
gion or too close to a bright star. In addition, we were 
able to identify 175 additional X-ray sources in the area 
of our mosaic not covered by the data of BAR08, i.e., we 
recover a total of 1023 out of 1325 X-ray sources from 
the catalog of Laird et al. For the the MIPS-24 /im and 
70 /im and radio surveys, we identify 10758, 868 and 590 
(out of 1122) sources, respectively. Noticeably, about 
20% and 54% of the sources detected in MIPS-24 /im 



and 70 /jm present a multiple IRAC counterparts (typi- 
cally 2-3) and —6% in X-ray and Radio. 

All of these surveys cover an area larger than the IRAC 
observations, but none of them fully cover the IRAC strip 
(although the X-ray and MIPS-24 /mi data cover >90% 
of the area) . Table [3] shows the total area covered by the 
X-ray, MIPS and radio surveys, the fraction that area 
in common with the IRAC mosaic, and the number of 
sources with an IRAC counterpart in our catalog. All 
the sources in the different surveys located within the 
area covered by IRAC are detected in our catalog (with 
the exception of a few objects too close to bright stars). 

Note that despite the larger area of the VLA obser- 
vations, this survey covers only ~40% of the region sur- 
veyed by IRAC. In fact the data are limited only to the 
upper region of the IRAC mosaic (5>53.10°; see Fig- 
ured]). 

5.3. Catalog reliability 

In this Section we analyze in more detail the sources 
with a poor spectral coverage to asses the reliability of 
the IRAC catalog as a function of the magnitude in the 
[3.6] band. 

First we test the reliability with the widely-used 
method of comparing the number of detections at faint 
levels in the original and a negative image. This test re- 
veals that only —1% of the sources up to [3.6]<23.75 mag 
are spurious. 

However, this procedure is conceived to detect faint 
spurious sources based on the assumption of a symmet- 
ric noise, whereas most of the spurious detections are 
concentrated around bright regions as a consequence of 
PSF and saturation artifacts and an excessive deblend- 
ing. Therefore, we chose to follow a different approach 
to analyze the reliability of the catalog based on the fact 
that faint [3.6, 4.5] detections lacking a measurable coun- 
terpart in other bands are potential candidates for spu- 
rious detections. 

As discussed in § [SJ we find a nearly isolated group of 
sources detected in N(band)<4. Most of these sources 
are IRAC-only detections (cyan histogram of Figure [6]) 
that, in the absence of low significance optical detections 
to validate them, could be artifacts of the source extrac- 
tion procedure. 

The visual inspection of galaxies with N(band)<4 re- 
veals a clear dichotomy. We find that these sources 
are either strongly clustered around bright stars and ex- 
tended sources (typically low-z galaxies), or uniformly 
distributed across the image. Most of the sources in the 
first group are easily identified as spurious, whereas for 
the isolated sources, their high SNR (—5-10) and the vi- 
sual inspection seem to favor that they are real. We 
dealt with the reliability of sources with N(band)<4 by 
considering spurious all detections in a 20" -30" radius 
region around the brightest objects ([3. 6] =16- 15 mag). 
Up to [3.6]<23.75, approximately half the N(band)<4 
sources are spurious by this criteria, accounting for less 
than a 2% of the total sample. Nevertheless, even for the 
more conservative scenario, the sum of all N(band)<4 
detections constitutes less than —4% of the sample in 
the main region. Finally, we also check the reliability of 
these sources by studying how many of them are simulta- 
neously detected at both 3.6 and 4.5 /im. Unfortunately, 
we find that some real detections are missed at 4.5 /im 
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due to the slightly worse image quality, while some of the 
spurious sources are detected in both images in regions 
where the density of artifacts is larger (i.e., around bright 
stars) . 

Summarizing, we conclude that the overall reliability 
of the catalog is very high (>97%) and that the con- 
tamination by spurious sources is strictly restricted to 
the surrounding areas of very bright sources. For the 
sake of completeness, we do not remove the sources with 
N(band)<4 from the sample. Instead, we include in the 
data catalog (Tabled both the number of bands in which 
the source was detected (N(band)) and a flag indicating 
if the source is located in the vicinity of a bright object 
(see 

5.4. Star-galaxy separation 

Following PG08, we used eight different photomet- 
ric and morphological criteria to identify stars in the 
merged photometric catalog: (1) the average of SEx- 
tractor star/galaxy separation parameter (stellarity) for 
all bands where the source is detected; (2) when avail- 
able, the stellarity parameter and FWHM of the source 
in the HST images is also considered as an independent, 
more reliable, criterion. (3,4) the IRAC-based color - 
magn itude criteria criter ia of lEisenhardt et al.l (|2004f ) 
and iBarmbv et al.l (|2004f ) . (5) the concentration pa- 
rameter (i.e., the difference of the [3.6] mag measured 
in a 3" radius aperture and the SExtractor MagAuto), 
(6) the color-col o r and (7) color-magnitude criteria of 
lEisenhardt et al.1 (|2004l this time b ased on opti c al an d 
NIR bands). (8) The BzK criteria of lDaddi et all (|2004h . 
In spite of using multiple stellarity criteria, the hetero- 
geneous band coverage makes difficult the simultaneous 
application of the 8 criteria (e.g., the HST data cover 
only ^40% of the full area). For this reason, we chose 
the bulk of the criteria to be based on IRAC or optical 
colors, which are available for the majority of the sources 
independently of the region in which they are located. 
We verify that at least 5 criteria can be estimated for 
86% of the sources with [3.6]<23.75. 

A source was identified as a star if 3 or more of the 
stellarity criteria are satisfied. Based on this method, 
we found 2913 stars. Among them, we are able to iden- 
tify 69% of the spectroscopically confirmed stars. The 
rest of them satisfy at least 1 or 2 of the stellarity crite- 
ria. Figure [5] illustrate s the efficiency of the IR AC mor- 
phological crite ria by lEisenhardt et al.l (|2004l ) and the 
BzK criteria of iDaddi et all (|2004T ). Red dots indicate 
sources with N(criteria)>3. Clearly, the BzK criterion 
is more efficient identifying stars over a wider magni- 
tude range. Note that, as discussed in PG08, an IRAC 
selected catalog down to [3.6]<23.75 only includes a mi- 
nor fraction of stars (less than 5% at all magnitudes), 
the majority of them at bright magnitudes (~40% at 
[3.6] < 18). At the faintest magnitudes, the overall dim- 
ming of the objects could lead us to identify some stars 
as galaxies due to the lack of applicable criteria. How- 
ever, we find that both the galactic and stellar number 
counts in the IRA C bands are in go od agreement with 
those presented in lFazio et al.l (|2004f) down to our limit- 
ing magnitude. From the relative distribution of galaxies 
and stars we also find that the fraction of stars makes up 
for only 3% of all sources at magnitudes fainter than the 
the median of the sample ([3.6]=22.4). 



6. DATA CATALOGS AND DATABASE ACCESS 

The catalog with the multi-band photometry for all 
the IRAC-3.6+4.5 /im selected sample in the ECS is pre- 
sented in Table [5j Additional information regarding the 
photometry of the sources and other properties described 
in the paper are given in Table [6] Furthermore, we also 
present here a web-based interface to access our database 
containing all the results presented in this and the com- 
panion paper. The web tool is publicly available for the 
entire astronomical community. 

As explained in § 13. 2\ the catalogs presented in this pa- 
per are restricted to the 76,936 sources with [3.6]<23.75, 
the ^85% completeness level, which count with accurate 
IRAC photometry (SNR>8). As explained in § El our 
sample is extracted from the area of the IRAC mosaic 
counting with a frame coverage larger than 20 (t exp ^4ks 
in the IRAC-3.6 /im band). Nevertheless, a deeper, al- 
though less complete, catalog without any magnitude re- 
striction can be accessed through the online database 
( ^6.3|) . In the following subsections we describe the con- 
tents of Tables [5] and [6] and we present Rainbow Naviga- 
tor. 

6.1. Tabled photometric catalog 
These are the fields included in Table [5] 

• Object: Unique object identifier starting with 
iracOOOOOl. Objects labeled with an underscore 
plus a number (e.g, irac000356_l) are those iden- 
tified as a single source in the IRAC catalog built 
with Sextractor, but deblended during the photo- 
metric measurement carried out with the Rainbow 
software (see § 14. ip . Note that, although the cata- 
log contains 76,936 elements, the identifiers do not 
follow the sequence iracOOOOOl to irac076185. This 
is because the catalog is extrated from a larger ref- 
erence set by imposing coordinate and magnitude 
constraints. The table is sorted according to this 
unique identifier. 

• a, 5: J2000.0 right ascension and declination in de- 
grees. 

• zspec: Spectroscopic redshift determination drawn 
from the DEE P 2 spectroscopic sur vey or the cata- 
log of LBGs of lSteidel et aIT(l2003h . 

• qflag: Spectroscopic redshift quality flag from 1 to 
4. Sources with qflag>3 have a redshift reliability 
larger than 80%. 

• FUV, NUV, it',...: Effective wavelengths (in 
nanometers), magnitudes and uncertainties (in 
the AB system) for each of the 30 photometric 
bands compiled for this paper. The band-passes 
are GALEX FUV and NUV, CFHTLS u*g'r'i'z', 
MMT u'giz, CFHT12k BRI, ACS V 606 and i 8U , 
Subaru R, NICMOS J 1W ,H 160 , MOIRCS K s , 
WIRC JK, CAHA- JK S IRAC 3.6-8.0, and MIPS- 
24 /im and 70 /im. The bands are sorted according 
to the effective wavelength of the filter. We refer 
to section §3 for details on the photometric mea- 
surement and error calculation. The magnitudes 
do not include zero-point corrections (see Paper 
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Fig. 8. — Example of two of the stellarity criteria used in this paper to separ ate galaxies from stars. Left: B-z versus z-K for all sources 
in the catalog. The dashed line depicts the BzK criteria of Daddi ct al. (2004) to isolate stars. The solid lines delimit the regions where 
high-redshift galaxies are located. The red dots indicate sources satisfying 3 or more stellarity criteria. Right: Comparison of the magnitude 
in a 3" radius aperture ([3.6 ]^//) and SExtrac tor MagAuto ([3.6]auto), i.e, the concentration parameter, versus [3.6]auto for all sources in 
the catalog. The criterion of Eiscnhardt ct al. (2004) identifies stars with -0.25<[3.6L//-[3.6]auto<0.2 and [3.6]auto<17.80. 



II, § 3.1.3). A value of the magnitude equal to - 
99.0 with error equal to 0.0 indicates a not valid 
photometric measurement. A negative error indi- 
cates that the source was undetected by Sextrac- 
tor but a positive flux was obtained when forcing 
the measurement at that position using the appro- 
priate aperture (see § 14. 2|) . An error equal to 0.0 
indicates that the source was undetected and the 
forced measurement returned a negative flux. In 
this case, the value of the magnitude is an upper 
limit equal to the sky-rms (lcr) in the photometric 
aperture (see § 14. 2ft . 

6.2. Tabled photometric properties catalog 
These are the fields included in Table [6] 

• Object: Unique object identifier (the same as in the 
photometric catalog in Tabic [5]). 

• a, 5: J2000.0 right ascension and declination in de- 
grees. 

• N (bands, detect): Number of UV-to-NIR in which 
the source is detected. 

• N (bands, forced): Number of UV-to-NIR in which 
the source is a priori undetected, but the forced 
photometry is able to recover a valid flux. 

• Flag: Quality flag indicating that the source is lo- 
cated in the vicinity of a bright object. Sources 
detected only in the IRAC bands (N(band)<5) 
and close to a bright ([3.6]>16) source are likely 
to be spurious (see § 15. 3[) . The values of the 
flag indicate: (5) source within 70"-100"of the 
brightess saturated stars in the field, (4) source 
within 30"of a [3.6] < 15 source, (3) source within 
20"of a 15<[3.6]<16 source, (2) source within 15"of 



a 16<[3.6]<17 source, (1) source within 10"of a 
17<[3.6]<18 source. (0) source un-flagged. 

• Stellarity: Total number of stellarity criteria satis- 
fied (see § 15. 4p . A source is classified as a star if it 
satisfies 3 or more criteria. (§ 15. 4p . 

• Region: Region of the field in which the source 
is located (§ [5j) : A value of 1 or indicates that 
the source is in the Main or Flanking region, re- 
spectively. The Main region is defined as the area 
of the IRAC mosaic within 52. 16<<5<53.20 and 
a>214.04, the Flanking regions are those contain- 
ing the remaining area. 

6.3. Rainbow Database & Navigator 

The photometric catalog presented here and the in- 
ferred stellar parameters discussed in Paper II are ob- 
tained using the Rainbow software (see PG05, PG08). 
The program includes different sub-routines for each 
task, and the output of each step serves as the input for 
the next. After the data processing, both the input (im- 
ages, spectra, templates) and the resulting catalogs are 
stored in a database with individual sources as building 
blocks. On doing so, we achieve several goals: (1) each 
object is fully characterized with all the available infor- 
mation; (2) the data can be easily sorted and retrieved 
according to several criteria; (3) the data are homoge- 
neously stored, which allows a straightforward combina- 
tion with data from other Rainbow fields and projects 
(such as those in PG05 or PG08). 

In order to provide worldwide access to the data stored 
in our database, we have developed a publicly available 
web interface, dubbed Rainbow Navigator^ . Here we 
briefly describe the main features of the utility. A de- 
tailed description of its capabilities can be found at the 
website. 

8 http://rainbowx.fis.ucm.es 
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Rainbow Navigator is essentially a user-friendly web 
interface to a database containing all the data products 
resulting from the process of creating and analyzing the 
IRAC-selected catalog presented in this two papers, from 
the initial source detection to the estimate of the stellar 
parameters. 

As many other astronomical query interfaces (e.g, 
NED, SIMBAD), Rainbow Navigator allows to retrieve 
the information for single sources searching for the source 
name or by coordinates. It also allows to create subsets of 
the complete catalog based on multiple constraints over 
the multi-band photometry, the redshifts or the stellar 
parameters. In addition, we have incorporated a cross- 
matching tool that allows to compare catalogs uploaded 
by the user to the IRAC-selected sample stored in the 
database, returning the sources in common. Rainbow 
Navigator also has an on-thc-fly utility to create sky 
maps of a selected area, including point-and-click access 
to the individual sources. 

Furthermore, an interesting feature that we do not in- 
clude in the public catalog for the sake of simplicity is 
the possibility of retrieving observed and rest-frame syn- 
thetic magnitudes over a predefined grid of 52 different 
filters covering the whole spectral range from the UV 
to the radio wavelengths. These values are computed 
by convolving the best fitting template (see Paper II) 
for each source with the appropriate filter transmission 
curve. 

Each source of the catalog has its own data sheets that 
provides all the available information including not only 
the photometry, the redshifts or the stellar parameters, 
but also detailed information of the stellarity, synthetic 
magnitudes, the multiplicity in other bands, jointly with 
the unique identification and coordinates of the counter- 
parts in each given band. The tool also provides postage 
stamps of the source in each of the available bands, that 
can be modified on-the-fly or combined to create simple 
false color images. Furthermore, each page includes a fig- 
ure depicting the full SED showing the fit of the optical 
and IR data jointly with the best fitting templates. 

Figures IWT21 show examples of the multi-band postage 
stamps, the UV-to-FIR SED and fitting templates, and 
the clickable map utility for a few galaxies at different 
redshifts. The postage stamps also illustrate the differ- 
ent photometric measurements in each band (aperture 
matched and circular apertures) and the forced detec- 
tions for very faint sources (see, e.g., the source in Fig- 
ured!]). 

Rainbow Navigator has been conceived to serve as a 
permanent repository for future versions of the catalogs 
containing improvements over the previous results (the 
present version is data release 1), and also to similar data 
products in other cosmological fields (such as GOODS-N 
and GOODS-S, presented in PG08). Currently, it pro- 
vides public access to the IRAC-selected catalog in the 
EGS presented in this paper, and also a similar release 
of the data described in PG08 for a small piece of the 
central region of the GOODS-S region. 

7. SUMMARY 

We presented an iRAC-3.6+4.5/mi selected 
catalog in the EGS characterized with multi- 
wavelength photometry. The sample contains 
76,936 sources with [3.6]<23.75 (85% complete- 



ness of the sample) covering an area of 0.48 deg 2 . 
The IRAC sources are characterized with FUV 
NUV u*g'r'i'z' u'gRiz BRI V 606 i 8 ii ^lio^ieo JK [3.6]- 
[8.0] photometry. In additio n, we have cross-c orrelated 
the sample with X-ray data (jLaird et al.ll2009L AEGIS- 
X), Spitzer/MIPS 24 /jm and 70 fim FIR photometry 
and VL A- 20cm radio data (AEGIS20; llvison et all 
2007). Secure spectroscopic redshifts are also included 
in the catalog for 7,636 sources with [3.6]<23. 7 5 ob- 
tained from the DEEP2 Survey and lSteidel et all (]2003l 
LBGs at z>3). The data described in this paper are 
publicly available, and will be part of future extended 
analysis and projects. The main results of this work are 
summarized below. 

• The extraction of the IRAC sample presented 
in this paper was limited to region with expo- 
sure times >4 ks. The average survey depth is 
texp ~10 ks. Aperture photometry was performed 
in the 4 IRAC bands simultaneously allowing us 
to obtain upper limit fluxes for undetected sources 
in the [5.8] and [8.0] bands. We removed spuri- 
ous detections masking areas around bright stars. 
The estimated 85% completeness level for point 
sources is [3.6,4.5]~23.75 and [4.5,5.8]~22.25. The 
3<T limiting magnitude estimated from the sky- rms 
is [3.6,4.5]~24.75 and [4.5,5.8]~22.90. We also val- 
idated the quality of the photometry by compar- 
ing our results with the catalog of Barmby et al. 
(2008), finding good agreement in both magnitudes 
(< 0.05 mag) and uncertainties f§ !3.3|) . Some small 
systematics were found in this comparison, which 
can be attributed to the slightly different reduc- 
tion versions and the limiting depth of the source 
extractions. 

• We described in detail our custom photometric pro- 
cedure Rainbow, developed to measure photometry 
in multi- wavelength data in a consistent way. The 
main steps followed by our photometry software 
are: (1) re-calibration of the local (5'x5') astro- 
metric solution for each pair of images, improving 
the accuracy of the cross-identification of sources 
to a limit below 0.1" (0.2") between optical-NIR 
(IRAC-ground based) images. (2) Deconvolution 
of blended IRAC sources (~16% of the entire cat- 
alog) separated by more than 1" , using the posi- 
tions of the optical/NIR counterparts as reference. 
The IRAC photometry of deblended sources is ac- 
curate to 0.03-0.10 mag, depending on the flux ra- 
tio between neighbors (§ 14. 1[) . (3) Measurement 
of consistent aperture matched photometry for a 
wide range of ground— and space— based observa- 
tions, with different depths and resolutions (§ 14. 2| ; 
(4) Obtaining (forced) flux measurements and up- 
per limits for faint, undetected (in a direct analysis 
of each image) sources (§ 14. 2[) . (5) Computing ro- 
bust photometric errors that take into account vari- 
ations in the sky rms and signal correlation (§ 14. 3|) . 

• The inhomogencous multi-wavelength coverage of 
the region covered by IRAC and the differences 
in depth of the various datasets justifies the split 
of the complete sample in two complementary re- 
gions. The main region (covering 0.35 deg 2 ), de- 
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limited by the footprint of the CFHTLS image, 
constitutes the region with the highest data qual- 
ity. Here the SEDs present a median coverage of 19 
bands, including high resolution imaging from HST 
(for —50% of the region) and deep NIR data from 
Subaru-MOIRCS (for -25% of the region). Nearly 
—70% of the complete sample presented in this pa- 
per is located in the main region. In the flanking 
regions (0.13 deg 2 ), the optical-to-NIR coverage is 
still robust, but the median band coverage is re- 
duced to 11 bands, lacking the high resolution HST 
data and the deepest NIR observations. 

• The overall detection efficiency of counterparts for 
the IRAC sources ([3.6]<23.75) in other bands is 
high: more than 85% of the sources are detected 
(with SNR>5) in the deepest optical data (R- or 
i'-bands) and 70% in the deepest (MOIRCS) K s - 
band images. Our method to perform (forced) 
photometric measurements for a priory undetected 
sources allows us to recover 10-20% additional 
sources in the shallowest images. Despite the 
large fraction of IRAC sources detected in all other 
bands, we find that ~2% of the sample is detected 
in IRAC only, some of them at relatively bright 
magnitudes 21<[3.6]<22 and with high SNR. 

• Around 10% of the sample counts with reliable 
spectroscopic redshifts. Nearly 20% and 2% of the 
sources are detected in the MIPS-24 /xm and MIPS- 
70 z/m data, respectively. This allows a detailed 
analysis of their IR-based SFRs. In addition, we 
recover 1023 of the X-ray sour ces in the catalog 
published by lLaird et all (|2009l 77% of their com- 
plete catalog) and 590 radio sources from the cat- 
alog published by llvison et all (|2007l 52% of their 
entire sample). 

• We presented the publicly available web-interface, 
Rainbow Navigator, to the database containing all 
the data products from this and the companion pa- 
per. The interface allows the access to the data 
using customizable queries on the photometry and 
derived parameters. Furthermore, it provides ad- 
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ditional capabilities to inspect the data, such as 
the creation of on-thc-fly clickable sky maps or the 
cross-match of the IRAC sample presented here to 
a user-provided catalog. We have made a signifi- 
cant effort to develop a useful and accessible tool 
that maximizes the legacy value of our catalogs in 
the EGS, and the future Rainbow based projects in 
other relevant fields (e.g., PG08). The online ver- 
sion of the catalog (IRAC-DR1) contains a deeper, 
although less complete, sample limited to [3.6] <25. 

• In Paper II, we make use of the very detailed UV- 
to-FIR SEDs presented here to estimate photomet- 
ric redshifts, stellar masses and SFRs. 
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Fig. 9. — Example of the results available for each IRAC source in the Rainbow database, accessible through the Rainbow Navigator 
interface. This source is irac089089, a galaxy at z=0.2170. Top-left, Map (2'x2') of the sky area around the source. The green dots in the 
map depict sources in our IRAC catalog with [3.6]<23.75. By clicking on any of the sources, the interface presents an individual webpage 
with all the available information for that source. Top-right, Full UV-to-FIR SED of the central galaxy in the sky map, a source at z=0.21. 
The UV-to-NIR data is fitted to a stellar population model, while the IR-part of the SED is fitted to the models of CE01, DH02. A summary 
of the estimated stellar parameters, such us the stellar mass or the global SFR is shown in the upper-left corner of the Figure. Middle 
panels: grey-scale postage stamps (with size 40"x40", except for the MIPS- 70 image, whose size is 2'x2', same as the map in the top-left 
figure) of the galaxy in some of the available bands, covering different wavelength ranges. The Kron aperture used to measure consistent 
photometry in optical/NIR bands and the circular aperture used in bands with significantly lower resolution are shown in all panels (red 
and blue, respectively). Bottom-left: RGB color stamp obtained by combining images in the ACS-V606: *814 bands. The Rainbow Navigator 
web interface allows to produce on-the-fiy monochromatic and RGB images changing the cuts interactively. Bottom-right: DEEP2 1-D 
spectra of the galaxy depicting some of the identified lines. The wavelength range and the redshift can be modified interactively. 




Fig. 10. — Same as Figure [9] for the galaxy iracl60152 at z=0.99. The band label shown in red in the GALEX-NUV postage stamps 
indicates a non-detection. For each of these bands cr^y is depicted in the UV-to-FIR SED as red arrow. Green name labels indicate forced 
detections (see § 14.21 1. These bands are shown as green stars in the SED plot. 
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Fig. 11. — Same as Figure [9] for iracl29479, an IR-bright galaxy at z=1.8. This galaxy correspond to source EGS11 in the paper of 
IHuane ct al. (2009). Band labels in red in the grey-scale postage stamps indicate non-detections. For each of these bands a^y is depicted 
in the UV-to-FIR SED as red arrow. Band labels in green indicate forced photometric measurements (see § 14.21 1. These bands are shown 
as green stars in the SED plot. 



25 



\est CM™] 

1 10 100 1000 




^ o 

1 ° 




10 100 1000 

\t,s [Mm] 




CFHTLS-i 








ACS-i 







10" 



10" 



IRAC-36 K 


IRAC-58 » 


IMIPS-24 


f rt 

♦ 1 0" 


} V 

1 

* 10" 


■r* 

10" 





Fig. 12. — Same as Figure [TT] for iracl64462, an LBG at z=2.8. Band labels in red in the grey-scale postage stamps indicate non- 
detections. For each of these bands cr^y is depicted in the UV-to-FIR SED as red arrow. Band labels in green indicate forced photometric 
measurements (see § 14. 21 1. These bands are shown as green stars in the SED plot. 
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The IRAC-3.6+4.5 /mm sample: Detection efficiency in each band and region. 



Band 



Aoff 



MAIN REGION (0.35 deg 2 ) 



FLANKING REGIONS (0.13 deg 2 ) 



Area 



[3.6]<23.75 [53,030] 



[3.6]<24.75 [77,607] 



Area 



[3.6]<23.75 [23,906] 



[3.6]<24.75 [35,416] 



stats 







\ a ) 




X-ray 1 " 


0.31,1.24 nm 


0.97 




FUV 


0.154 


0.99 


27 


NUV 


0.232 


0.99 


26 


u' 


0.363 


1.00 


25 


u* 


0.381 


1.00 


25 


CFH-B 


0.440 


1.00 


25 


g 


0.481 


1.00 


25 


9' 


0.486 


1.00 


24 


V606 


0.592 


0.56 


24 


r' 


0.626 


1.00 


24 


R 


0.652 


1.00 


24 


CFH-R 


0.660 


1.00 


24 


i' 


0.769 


1.00 


24 


i 


0.782 


1.00 


2d 


*S14 


0.807 


0.56 


24 


CFH-I 


0.813 


1.00 


24 


z' 


0.887 


1.00 


2.3 


z 


0.907 


1.00 


23 


JllO 


1.10 


0.04 


23 


Q2k-J 


1.21 


0.40 


23 


WIRC-J 


1.24 


0.42 


22 


#160 


1.59 


0.04 


24 


n' - k s 


2.11 


0.46 


22 


Ks 


2.15 


0.26 


24 


wirc-k 


2.17 


0.97 


21 


IRAC-36 


3.56 


1.00 


22 


IRAC-45 


4.51 


1.00 


22 


IRAC-58 


5.69 


1.00 


22 


IRAC-80 


7.96 


1.00 


22 


MIPS-24 


23.84 


0.91 


19 


MIPS-70 


72.49 


0.80 


15 


VLAt 


2E06 


0.74 





(4) 



ALL/SNR>5(%) 
(5) 



stats 



(6) 



ALL/SNR>5(%) 
(7) 



stats 



(8) 


0.93 




0.15 


27 


0.15 


26 


1.00 


25 


1.00 


24 


1.00 


25 



(9) 



ALL/SNR>5(%) 
(10) 



stats 



(11) 



ALL/SNR>5(%) 
(12) 



1 



717 (1.3%) 
58.8/7.9 
74.2/25.8 
88.4/52.8 
88.2/53.9 
89.6/55.2 
96.1/85.1 
95.4/83.6 
94.7/83.5 
97.1/89.3 
97.0/88.5 
95.3/71.5 
95.3/85.4 
97.4/90.3 
96.7/80.7 
94.1/68.6 
93.3/66.2 
96.0/77.7 
84.6/65.4 
93.7/38.5 
62.8/16.6 
92.3/80.8 
75.8/6.7 
95.5/73.8 
62.6/27.1 

100.0/99.2 
99.8/95.3 
89.7/41.1 
85.8/34.4 
31.5/20.2 
13.1/1.5 

380 (<1%) 



28.0321 
26.60|| 
25.85|| 
25.60|| 
25.33?! 
25.44 2 | 
25.20|| 

24.92;! 
24.851 



24.801 
24.85|| 
24.41 2 | 
24.511 

24.52^ 
24.45|§ 
24.231 



24.44; 
24.25! 
23.74! 
22.49! 
25.18; 
22.74! 
25.14! 
21.82; 
23.18! 
23.26! 
22.89! 
22.97! 
19.7L 
15.81 



718 (1.3%) 
57.4/6.1 
70.8/20.2 
86.1/44.1 
85.2/44.2 
87.1/45.6 
95.0/80.3 
94.0/78.2 
93.9/78.7 
96.0/84.2 
95.8/82.5 
93.1/60.2 
94.2/78.3 
96.1/83.8 
94.7/69.3 
90.1/54.5 
89.1/51.2 
92.7/64.9 
88.9/55.6 
88.4/27.1 
46.1/11.7 
91.7/69.4 
67.0/4.6 
89.8/52.9 
42.6/16.3 

100.0/84.5 
98.9/75.3 
78.4/28.0 
73.8/23.5 
29.7/18.1 
7.5/5.7 

380 (<1%) 



0.80 
1.00 
1.00 
1.00 
1.00 
0.93 
0.86 
0.42 



a 23.76 



00 
00 


24 
24 


3611 
36i 


00 


24 


13i 


00 


23 


971 


00 


23 


95|2 


32 


22 


01 23 
ol 20 


26 


22 


17i 



91 4922.50 
19-86 
22.52|;« 
99 fio23!6l 

22.45sHjS 
22.51S? : l2 



19.36 
15.57 



im 



307 (1.2%) 



7 


27.83 


.7 


26.57 


.3 


25.69 


.1 


25.27 


.9 


25.49 



45.3/18.4 
100.0/97.8 
99.3/91.1 
87.1/38.6 
82.6/32.9 
29.9/19.1 
13.6/1.2 
210 (<1%) 



mi 
mi 



24.14 



21.45 
23.17! 



23. 27? 



22. 80? 



22.83 
19.49- 
15.57 



307 (1.2%) 
58.3/6.8 
72.5/23.7 
88.0/57.8 

83.7/43.0 
95.5/80.5 



97.4/87.7 
91.6/73.7 


24 ' 74 23'27 

24-7611 


96.5/81.9 
88.4/63.0 


96.4/81.1 


24.541:17, 


94.6/70.2 


89.2/61.1 


24.43i : 70 


84.1/47.5 


97.5/85.5 


24.44i: 6 l 


95.3/75.4 


57.9/20.4 


99 qc23.46 
Zz - Od 20.73 


41.1/14.1 


67.3/7.3 


99 on23.76 


59.6/5.0 



42.4/14.4 
100.0/80.6 
97.1/70.5 
75.5/26.4 
71.1/22.6 
25.2/15.1 
6.7/4.7 
210(<1%) 



[3.6]<23.75 [53,030] [3.6]<24.75 [77,607] [3.6]<23.75 [23,906] [3.6]<24.75 [35,416] 



Rcdshiftst 0.64-0.91 0.94 6191 12.0% 6420 8.4% 0.84 1445 6.0% 1481 4.2% 

NOTE. — Detection efficiency for the sources IRAC-3.6+4.5/im sample in the different bands compiled for this work. The analysis is divided in 2 zones: the main region, defined as the overlapping area 
between the IRAC and the CFHTLS mosaics (52.16° < 6 <53.20° and a >214.04°), and the flanking regions. 

| For these catalogs, we quote only the number of sources detected in IRAC, and the percentage of total IRAC sample that they represent. 

(1) Band name. 

(2) Effective wavelength of the bands in microns. 

(3) Fraction of the main region covered by the observations in each band. 

(4.6) Median and quartiles of the magnitude distribution in each band for main region sources with [3.6]<23.75 and [3.6]<24.75, respectively. 
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(8) Same as (3) in the Flanking regions. 
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The IRAC-3.6+4.5 /jm sample: Multi-band photometry 
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5689.4 
22.385 
0.403 



781.5 
22.862 
0.040 
7957.6 
22.378 
0.189 
781.5 
17.870 
0.029 
7957.6 
19.341 
0.046 
781.5 
22.463 
0.036 
7957.6 
21.247 
0.067 
781.5 
19.985 
0.029 
7957.6 
19.695 
0.045 
781.5 
23.433 
0.048 
7957.6 
22.479 
0.383 



807.3 

-99.000 
0.000 

23844.0 
19.545 
0.238 
807.3 

-99.000 
0.000 

23844.0 

-99.000 
0.000 
807.3 

-99.000 
0.000 

23844.0 
18.860 
0.118 
807.3 

-99.000 
0.000 

23844.0 
19.436 
0.165 
807.3 

-99.000 
0.000 

23844.0 

-99.000 
0.000 



813.3 
22.824 
0.034 
72493.7 
-99.000 
0.000 
813.3 
17.928 
0.029 
72493.7 
-99.000 
0.000 
813.3 
22.390 
0.032 
72493.7 
-99.000 
0.000 
813.3 
19.978 
0.030 
72493.7 
-99.000 
0.000 
813.3 
23.539 
0.067 
72493.7 
16.266 
1.282 



NOTE. — (1) Object unique identifier in the catalog. The catalog is sorted by this field. 
(2,3) Right Ascension and Declination (J2000) in deg rees. 

(4) Spectroscopic redshift determination drawn from i Davis ct al. 2007; ~8,000 galaxies) and i Steidcl et al. 2003; LBGs at z>3). 

(5) Quality flag of the spectroscopic redshift (4— >99.5%, 3 = >90%, 2 = uncertain, l = bad quality). Only redshifts with qflag>2 have been used in the analysis. 

(6-50) Effective wavelength of the FUV NUV u* g' V ' i' z' u' gB.iz BRI 1^606*814 ^110^160 J K [3-6]-[8.0] [24] [70] bands in nanometer. Observed magnitude with the associated uncertainty as measured in the apcrtured matched 
procedure (sec § \4\ . All magnitudes refer to the AB photometric system. 

Negative values of the photometric uncertainty are for forced photometric detections. Magnitudes >0 with uncertainty equal to 0.00 are non-detections; for these bands the photometry indicates the magnitude of cr^y (see jj |4.2| - 
(This table is available in its entirety in a machine-readable in the online version. A portion is shown here for guidance.) 



to 
-J 



28 



TABLE 6 

The IRAC-3.6+4.5 fj,M sample: Photometric properties 



Object 


a 


5 


N(bands, detect) N(bands, forced) Flag 


Stellarity Region 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


irac003270.1 


215.43910540 


53.08468920 


18 


1 





1 


irac003278 


215.42614011 


53.09447161 


18 


2 





7 1 


irac003291_l 


215.44058360 


53.08123980 


18 


2 





1 


irac003310 


215.42129738 


53.09430607 


19 


1 


1 


1 


irac003313 


215.43553774 


53.08200958 


17 


1 





1 1 



NOTE. — (1) Object unique identifier in the catalog. The catalog is sorted by this field. 
(2,3) Right Ascension and Declination (,12000) in degrees. 

(4) Number of optical-to-NIR bands (with effective wavelengths below 8.0"m) in which the object is detected. 

(5) Number of optical-to-NIR bands (with effective wavelengths below 8.0/^m) in which the object is a priori un-dctcctcd. but the forced photometry (sec ^ |4.2| recovers 
a valid flux. 

(6) Quality flag indicating the proximity of a very bright source in the vecinity of the source. Sources detected in very few bands (N(band)<5) and located nearby a 
bright source are likely to be spurious detections (see|[37SJ. 

(7) Sum of all the stellarity criteria satisfied (see § [5*4] . A source is classified as star for Stellarity > 2 . 

(8) Region of the field in which the source is located: 1 for the Main region, for the Flanking regions. 

(This table is available in its entirety in a machine-readable in the online version. A portion is shown here for guidance.) 



